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ABSTRACT 


Extensive  results  of  numerical  calculations  for  turbulent 
recirculating  flew  with  heat  transfer  are  presented.  The 
geometries  studied  were  two-dimensional  flows  over  rearward 
(backward)  -  facing  steps  and  sudden  pipe  expansions.  The 
confutations  described  were  performed  using  the  k-e,  two  equation 
model  of  turbulence.  TO  yield  better  predictions  of  heat  transfer, 
the  standard  model  was  improved,  by  developing  algebraic 
approx ilnations  for  Reynolds  stresses  and  scalar  fluxes,  so  as  to 
account  for  the  effects  of  streamline  curvature,  pressure-strain 
and  variation  of  turbulent  Prandtl  number  in  the  flow  field.  In 
addition,  the  low  'Reynolds-number  version  of  the  standard  k-e  model 
was  employed  in  the  redeveloping  region.  In  order  to  apply  the 
low-Reynolds-nuraber  model  to  the  redeveloping  region,  a  two  pass 
procedure  was  developed  to  partition  the  flow  field  into  separated 
and  reattached  regions.  This  procedure  consists  of:  a  first  pass 
utilizing  the  standard  k-e  model  to  establish  the  reattachment 
length,  and  a  second  caiputational  pass,  in  which  all  variables  are 
recomputed,  rxjring  the  second  pass,  the  standard  k-E  model  is 
maintained  up  to  reattachment  but  the  low-Reynolds-minber  version 
is  employed  beyond  reattachment.  The  numerical  scheme  incorporates 
the  basic  solution  algorithm  of  the  TEACH-T  code.  However,  this 
algorithm  was  extended  to  include  heat  transfer  and  was  improved  by 
using  quadratic  differencing  for  the  convective  terms. 
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Results  of  calculations  for  mean  velocity,  mean  temperature 
and,  turbulent  intensity  fields,  together  with  heat  transfer 
coefficients  for  both  the  side  and  base  walls  are  presented.  Also 
described  are  the  effects  of  the  turbulent  approach  boundary  layer 
on  heat  transfer,  the  characteristics  of  the  redeveloping  turbulent 
thermal  boundary  layer  after  reattachment  and  correlations  of  heat 
transfer  coefficients. 

Comparisons  are  made  with  other  numerical  calculations  and 
with  available  experimental  data.  The  results  indicate  that 
reliable  predictions  of  reattachmenfc  lengths  can  be  obtained  by  the 
present  improved  standard  k-e  model.  Also,  fairly  accurate  overall 
heat  transfer  coefficients  are  obtained  by  incorporating  both  the 
unproved  standard  fc-t  model  anu  the  low-BeynoluS  number  form  of  the 
model  in  a  two-pass  procedure. 
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CHAPTER  1 


INTRODUCTION 

1.1  General  Remarks 

In  recent  years,  the  energy  crisis  has  motivated  heat  transfer 
engineers  to  develop  high  performance  heat  exchangers  and  other 
advanced  thermal  devices.  In  developing  the  new  technology  needed 
to  acconplish  this  task,  a  great  deal  of  attention  has  been  given 
to  the  various  heat  transfer  enhancement  techniques.  Some  of  these 
augmentation  schemes  include:  surface  promoters ,  vortex  flews,  and 
additives  [1] .  With  surface  promoters,  the  heat  transfer  for  flows 
over  rough  surfaces  can  be  as  much  as  15  percent  larger  than  for 
flows  over  aiooth  surfaces  12,-3]*  In  vortex  flows,  the  heat 
transfer  is  significantly  increased  (by  swirling  the  flow),  in  some 
cases  as  much  as  80  percent  14,5].  Ihe  present  research  is  focused 
on  effects  due  to  sudden  changes  of  the  surface  geometry.  Such 
abrupt  changes  can  enhance  the  heat  transfer  by  as  much  as  50 
percent  over  that  of  the  corresponding  attached  flow  [6] . 

The  heat  transfer  augmentation  which  occurs  as  a  result  of 
alterations  of  the  flow  caused  by  sudden  changes  in  body  geometry 
is  a  very  important  aspect  of  engineering  heat  transfer.  Such  flow 
alteration,  consisting  of  separation,  reattachment  and  subsequent 
redevelopment  may  accidentally  occur  in  seme  systems,  and  in 
others,  may  be  deliberately  induced,  primarily  to  enhance  the  wall 
convective  heat  transfer  rates  at  the  expense  of  an  increase  in  the 
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friction  factor.  While  ti.ere  have  been  significant  contributions 
in  this  field,  current  understanding  of  the  reattachment  process  is 
still  relatively  poor,  essentially  because  of  its  complexity. 
Thus,  a  designer  is  still  somewhat  at  a  loss  in  predicting  heat 
transfer  in  separated  flows.  The  geometries  chosen  for  the  present 
calculations  are  the  sinpie,  yet  vitally  important,  two-dimensional 
back-step  and  sudden  expansion  in  a  pipe.  A  considerable  amount  of 
the  needed  information  for  validating  turbulence  models  and 
numerical  schemes  can  be  obtained  by  studying  these  fundamental 
engineering  geometries. 

Among  two-dimensional  flows,  the  back-step  is  the  simplest 
reattaching  flow  of  practical  significance.  The  separation  point 
is  fixed  at  the  edge  of  the  step,  and  there  is  only  one  separated 
zone,  unlike  flow  over  a  fence  or  obstacle.  It  should  be  noted 
that,  although  the  overall  gross  features  of  the  turbulent  flow  are 
two-dimensional,  the  fluctuating  components  are  actually 
three-dimensional.  Even  in  this  "simplest"  reattaching  flow,  the 
flow  field  for  the  back-step  is  still  very  complex. 

The  qualitative  nature  of  separated  flew  fields  has  been  knewn 
for  a  long  time  as  discussed  by  Prandtl  and  Tietjens  17] .  Figure 
(1)  shews  a  typical  abrupt  expansion  geometry.  The  inception  of 
separation  occurs  at  the  point  of  geometrical  discontinuity.  The 
’dividing  streamline'  is  the  hypothetical  streamline  that  separates 
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the  reversing  fluid  from  that  continuing  in  the  main  flow 
direction.  Eeattachment  is  the  point  where  the  separated  layer 
rejoins  a  boundary  surface  (or  another  layer  as  in  free  shear 
flows) .  Thus,  at  reatcachment  two  flows  of  opposing  directions 
exist;  and  the  shear  stress  vanishes.  Figure  (1)  depicts  sene  of 
the  qualitative  flew  characteristics.  Immediately  after  the 
boundary  layer  separates  at  the  step,  a  new  free-shear  layer 
develops  inside  the  original  shear  layer.  In  the  separated  zone 
the  flow  is  highly  turbulent  and  is  reversed  towards  the  step  in 
the  near-wall  region  and  re-enters  into  the  new  shear  layer.  This 
results  in  a  recirculating  zone  that  is  separated  from  the  main 
stream  flow.  Also,  in  the  separated  region  the  shear  layer  curves 
sharply  downwards  and  is  influenced  by  the  effects  of  curvature, 
adverse  pressure  gradient,  and  wall  interaction.  Immediately 
downstream  of  reattachment,  a  new  sub-boundary  layer  is  formed; 
this  layer  interacts  with  the  original  (incoming)  shear  layer. 
Further  downstream,  the  flew  tends  to  a  more  reestablished 
turbulent  boundary  layer  (or  fully-developed  flow  in  a  pipe  or 
channel)  [8,9]. 

In  the  recirculation  region,  the  early  studies  of  Abbott  and 
Kline  [10]  clearly  identify  the  presence  of  three  dissimilar  zones. 
Iirenediately  downstream  of  the  step  is  a  three-dimensional  zone, 
characterized  by  one  or  more  vortices  rotating  about  an  axis  normal 
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to  the  plane  of  the  flow.  Adjacent  vortices  are  counter-rotating, 
and  they  may  not  all  be  of  the  same  size.  Downstream  of  the  first 
zone  is  a  two-dimensional  zone  that  contains  the  classical 
recirculation  pattern,  i.e.,  flew  near  the  wall  moving  upstream  and 
near  the  dividing  streamline  moving  downstream.  Beyond  the  second 
zone  exists  a  transient  tail  region  (close  to  the  wall  and  upstream 
of  reattachment.)  that  eventually  disappears  with  the  second  zone. 

The  separation  mechanism  is  somewhat  different  for  a  body  with 
a  sudden  geometric  change  as  epposed  to  one  with  a  gradual 
geometric  change.  For  a  body  with  a  slowly  varying  geometric 
change  the  curvature  of  the  velocity  profile  at  the  wall  is 


directly  affected 

by  tiie 

xi  tijj  w  e  S 

pressure 

gradient. 

For 

instance,  for  a 

negative 

(favorable) 

pressure 

gradient. 

the 

velocity  profile  has  a  negative  curvature  all  the  way  to  the  edge 
of  the  layer,  and  no  separation  occurs.  Fbr  a  zero  pressure 
gradient,  the  velocity  profile  has  zero  curvature  at  the  wall  and 
negative  curvature  near  the  edge  of  the  layer,  and  no  separation 
occurs.  For  a  positive  (unfavorable  or  adverse)  pressure  gradient, 
the  velocity  profile  has  a  positive  curvature  at  the  wall  and 
negative  curvature  at  the  out  r  edge,  and  separation  occurs,  Fbr  a 
sudden  expansion  geometry  due  to  the  interaction  there  is  no 
impressed  pressure  gradient  as  such,  and  the  wall  (surface)  static 
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pressure,  expressed  in  terms  of  pressure  coefficients,  is 
characterized  as  follows.  At  separation  the  pressure  coefficient 
is  negative  and  decreases  gradually  to  a  minimum  value  just  beyond 
the  step.  In  the  reattachment  zone  it  increases  sharply  to  a 
positive  value,  reaching  a  maximum  just  beyond  reattachment. 
Further  downstream  of  reattachment  it  decreases  gradually  to  the 
value  of  the  local  free  stream  pressure  [9] .  'The  static  pressure 
variation  across  the  expansion  (perpendicular  to  the  flow 
direction)  is  nearly  constant  throughout,  except  for  sharp  peaks 
across  the  dividing  streamline. 

Hie  presence  of  a  recirculating  wall  layer  greatly  influences 
the  process  of  energy  transfer.  It  has  been  demonstrated 
experimentally,  cf  [6,8],  that  within  the  recirculating  region  the 
heat  transfer  coefficient  is  many  times  larger  than  that  of  an 
attached  flow.  This  increase  in  heat  transfer,  in  principle,  can 
be  attributed  to  the  increase  in  streamwise  turbulent  kinetic 
energy  in  the  mixing  layer  (region  remote  from  the  wall).  Recent 
measurements  of  Moss  and  Baker  I9j  and  Fraser  and  Siuuiy  [11] 
illustrate  this  effect.  They  showed  that  the  streamwise  turbulent 
intensity  (u  VUo")  is  small  within  the  buffer  zone  and  then 
increases  quite  sharply  to  a  maximum  in  the  mixing  layer,  before 
falling  off  again  to  a  low  value  in  the  main  flow.  Furthermore, 
the  peak  value  of  the  turbulent  kinetic  energy  increases  initially 
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with  the  downstream  distance  from  the  sudden  change  in  body 
geometry,  reaching  a  maximum  just  before  the  reattachment  point, 
and  then  decreasing  after  attachment.  In  the  redevelopment  region 
where  the  flow  tends  toward  a  reestablished  boundary  layer  (or 
fully  developed  flow  in  a  pipe  or  channel),  the  turbulent  kinetic 
energy  takes  on  boundary  layer  type  profiles.  The  high  values  of 
the  turbulent  kinetic  energy  in  the  mixing  layer,  particularly 
around  reattachment,  is  associated  with  spreading  of  the  layer  or 
an  increase  in  turbulent  diffusion  coefficient  in  the  main  flow. 
This  then  results  in  a  'shrinking'  of  the  near-wall  viscous 
sub-layer,  thus  allowing  greater  passage  of  heat  by  molecular 
diffusion.  The  increase  in  heat  transfer  due  to  recirculation  can 
also  be  viewed  as  a  manifestation  of  streamline  curvature.  For 
instance,  the  fluid  particle  near  the  wall  follows  a  longer  curved 
streamline  (as  opposed  to  the  rectilinear  path  present  in.  attached 
flows) ,  thereby  transferring  more  energy  along  the  curved  path. 

From  the  above  discussion  it  is  evident  that  the  separating, 
reattaching,  and  redeveloping  shear  layers  are  characterized  by 
active  coupling  between  the  external  flow  and  the  separated  shear 
layers.  It  is  for  this  reason  that  the  flew  field  is  complex .  in 
attached  flews,  active  coupling  between  the  free  stream  and 
boundary  layer  need  not  be  considered,  and  so  inviscid  and  boundary'- 
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layer  theories  can  usually  be  applied  in  an  uncoupled  manner  to 
yield  fairly  reasonable  predictions.  With  separation,  boundary 
layer  approximations  are  not  valid,  and  with  the  inclusion  of  heat 
transfer,  the  problem  of  coupling  becomes  even  more  carplex. 
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1.2  General  Form  of  the  Governing  Equations 

It  is  accepted  that  the  Navier-Stokes  equations,  comprising  a 
closed  set  of  equations  for  instantaneous  velocity  and  pressure 
fields,  are  valid  for  both  laminar  and  turbulent  flows.  For 

laminar  flows,  with  well  defined  initial  and  boundary  conditions, 
numerical  solutions  can  be  obtained  fairly  easily.  For  turbulent 
flows,  the  important  components  of  the  turbulence  pheixmena  take 
place  in  small  eddies  of  the  order  of  1  mm.  in  size  while  the 
entire  flow  field  may  extend  over  meters.  Reliable  numerical 

results  can  only  be  expected  if  the  numerical  discretization  is 
smaller  than  the  smallest  scales  of  turbulence.  Such  a  procedure 
would  require  a  large  and  impractical  number  of  nodes;  and  as  such, 
demand  computer  storage  and  time  far  beyond  the  capabilities  of 
present  systems.  As  an  example,  an  eddy  size  of  1  nan.  and  a  value 
of  10  percent  of  the  average  velocity  for  the  turbulent 
fluctuations,  will  result  in  turbulence  frequencies  of  10~5  sec  -1. 
Gebeeci  and  Smith  [12]  explored  the  possibility  of  obtaining  a 

computer  solution  to  reveal  the  fine  scale  structures.  Uiey  showed 
that  boundary  layer  flow  calculations  over  an  airfoil,  with  about 
ten  grid  points  per  eddy-size,  will  need  about  35,000  years  of 
computing  time  with  a  present-day  computer.  By  requiring  less 
physical  detail,  so  as  to  obtain  solutions  of  practical  engineering 
problems,  restrictions  can  be  made  regarding  the  fine  scale 
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structures.  The  problem  can  be  effectively  side-stepped  by 
considering  the  time-averaged  equations  of  turbulence.  Time  - 
averaging  'smears  out'  tho  fine  scale  structures,  and  the  resulting 
gross  properties  of  turbulence  (such  as  mean  velocities  and 
tenperatures)  are  much  more  gradual  in  spice;  as  a  consequence 
excessively  fine  grids  are  not  required.  However,  time-averaging 
the  governing  equations  results  in  terms  involving  statistical 
correlations  of  fluctuations  in  velocities,  temperature,  pressure, 
etc.  There  is  no  direct  way  of  knowing  these  correlations,  so  they 
have  to  be  modeled  in  terms  of  quantities  that  can  be  determined. 
Hence,  turbulent  modeling  involves  developing  a  set  of  equations 
for  the  correlations,  which  when  solved  with  the  mean  flow 
equatioas,  simulate  the  behavior  of  the  actual  flow.  The  details 
of  modeling  procedures  are  discussed  in  a  later  chapter. 

Rather  than  using  the  time-averaging  procedure,  there  are  less 
restrictive  formulations  in  the  literature.  A  group  of 
theoreticians,  such  as  Kraichnan,  Edwards  and  Qrszag  (these  studies 
are  discussed  in  the  bock  of  Orszag  [13] ) ,  have  developed 
sophisticated  statistical  theories  of  turbulence.  The  mathematical 
formulation  takes  into  account  the  dynamics  of  the  turbulent  motion 
through  the  interaction  of  eddies  of  differ?'.'  sizes.  The 
procedure  is  often  referred  to  as  the  "direct  interaction"  method 
and  is  discussed  by  Orszag  i!3J .  The  mathematical  ecnplexity  of 
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the  procedure  is  overwhelming;  thus  the  understanding  and 
applications  to  basic  engineering  calculations  of  turbulent  flow 
are  almost  non-existent.  Another  scheme  is  the  "subgrid-scale" 
procedure  of  Deardorff  [14]  in  which  averaging  is  taken  over  small 
time  intervals.  In  the  final  procedure  all  the  energy-containing 
turbulent  motions  are  treated  as  though  they  are  unsteady 
three-dimensional  flows.  However,  applying  this  procedure  to  study 
atmospheric  turbulence,  Deardorff  115]  reported  central  processor 
time  on  a  CDC  7600  of  about  350  hours;  thereby,  implying  the 
irapracticality  of  the  procedure  for  rudimentary  engineering 
calculations. 

-  in  the  present  work  the  generalized  steady  equations  of  fluid 
motion  in  the  Cartesian  tensor  notation  are  employed.  Details  of 
the  foregoing  analysis  are  given  in  Hinze  [16] . 

For  an  incompressible  fluid,  the  continuity  equation  takes  the 

form, 


3ui 

3^  "  0 


where  is  the  i'th  component  of  velocity.  Pur  a  Newtonian  fluid 
of  uniform  density,  p  ,  and  kinematic  viscosity,  v  ,  the 
instantaneous  Navier-Stckes  equations  in  the  absence  of  body  forces 
are, 


i 

P  3xj_ 


Hi  ) 

3Xj  / 


+  v 


(1.2.2) 
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Hie  Navier-Stokes  equations  for  turbulent  flow  is  obtained  by 
employing  the  Reynolds  procedure,  namely 

ui  **  ui  4  ui* 

P  *  P  +  p’ 

where  the  primed  and  overbar  superscripts  indicate  fluctuating  and 
mean  quantities,  respectively.  Before  substituting  the  Reynolds 
procedure  into  the  above  equations  to  obtain  the  Navier-Stokes 
equations  for  turbulent  flow,  a  brief  discussion  of  time-averaging 
with  respect  to  "steady"  turbulence  is  presented. 

A  turbulent  flow  is  referred  tb  as  steady,  if  for  all 
dependent  variable,  $ (x^,t) , 

c0 

$  (*£>  -  i  J*  $  (xi>t)  dt 
c0  0 

as  time  becomes  large  compared  with  the  time  scales  of 
turbulence  motions,  <p  ' .  Hie  average  value  is  independent,  of  the 
origin  of  t  of  the  averaging  procedure,  provided  t<tQ. 

Substituting  the  above  Reynolds  relations  into  Bgn.  (1. 2.2), 
making  use  of  Eqn.  (1.2.1),  and  time  -  averaging  the  results  yield 
the  corresponding  Navier-Stokes  equations  for  turbulent  flow, 
p  uj  -  -  li-  +  - —  I  p  -  P  1  (1,2.3) 

4  aXj  3xt  hxj  >■  3Xj  4  1 

The  equation  governing  the  transport  of  thermal  energy 


(Temperature,  T),  is. 
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u.  !L  s  S  +  r  - —  (  —  )  (1.2.4) 

J  3Xj  3Xj  3Xj 


where  Sp  is  the  generation  of  T  per  unit  volume  and  r  is  the 
kinematic  diffusivity  of  T. 

Using  the  same  Reynolds  procedure,  i.e., 

t  *  t  +  e' 


where  8'  is  the  fluctuating  part  of  the  tenperature.  Substituting 
into  Eqn.  (1.2.4)  and  time  averaging,  results  in. 


u,^_=<L  +  r3 

3  3x. 


(  3T  )  8 
9x-i  9x-i 


(Uj6) 


(1.2.5.) 


■3  3xj 

The  turbulent  results,  Dgns.  (1.2.3)  and  (1.2.5)  are  identical  to 
their  laminar  counterpart,  Eqns.  (1.2.2)  and  (1.2.5),  except  for 
the  unknown  correlations,  uj_uj  and  u^e' .  These  are  commonly 
referred  to  as  Reynolds  stresses  and  scalar  fluxes,  respectively. 


Solving  tor  these  unknown  stresses  and  fluxes  is  referred  to  as 


turbulence  closure  or  turbulent  modeling.  For  general  information 
.Appendix  (A)  discusses  the  mathematics  and  physics  involved  in  seme 
of  the  most  ectanon  turbulence  closure  techniques. 

Presently,  it  is  possible  to  predict  velocity  profiles  in  an 
inconpressible  turbulent  boundary  layer  with  a  Treasonable  degree  of 
success  by  using  a  wide  variety  of  eddy  viscosity  modsls  [12,17] . 
In  these  models  the  turbulent  stresses  in  the  mean  momentum 
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equation  are  written  in  terms  of  a  turbulent  (eddy)  viscosity  and 
mean  velocity  gradients.  This  level  of  success  has  led  many 
investigators  to  consider  a  similar  approach  for  the  prediction  of 
heat  transfer  in  thermal  boundary  layers  by  employing  the  Reynolds 
analogy  118,19].  The  hypothesis  is  that  the  same  complex  processes 
govern  the  transport  of  both  momentum  and  thermal  energy  in  a 
turbulent  boundary  layer.  It  follows  then  that  the  turbulent  heat 
flux  (e.g.  vV)  in  the  above  thermal  energy  equation  can  be 
written  as  an  eddy  conductivity  times  mean  tenperature  gradients 
with  the  turbulent  Prandtl  number  defined  as  the  ratio  of  turbulent 
viscosity  to  turbulent  conductivity.  The  turbulent  viscosity  can 
be  approximated  by  using  a  single  length  scale  (Prandtl* s  mixing 
length)  or  by  using  length  and  velocity  scales  (k-e  model).  Then, 
according  to  the  Reynolds  analogy,  the  turbulent  Prandtl  number  is 
0(1). 

The  focus  of  the  present  analysis  is  on  numerical  calculation 
of  the  flow  and  tenperature  fields,  turbulent  variables,  and  side 
and  bottom  wall  heat  transfer  coefficients,  for  both 
rearward-facing  steps  and  pipe  expansions,  for  the  entire  turbulent 
flow  regime.  The  range  of  Reynolds  numbers  emphasized  are  in  the 
lew  to  moderate  range.  The  upstream  wall  is  heated  and  the 
momentum  and  thermal  boundary  layers  are  turbulent  at  both 
separation  and  reattachment  (often  referred  to  as 
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turbulent-turbulent  interaction) .  Emphasis  in  this  study  is  on 
exploring  turbulent  viscosity  models  for  heat  transfer  flews  with 
recirculation  and  subsequent  redevelopment.  In  the  present  closure 
scheme  the  turbulent  viscosity  and  thermal  diffusivity  are 
approximated  by  less  restrictive  formulations  than  are  usually 
employed.  Generalized  functional  relationships  for  critical 
parameters  are  algebraically  derived  to  better  represent  the 
physics  in  both  the  recirculating  and  redeveloping  regions.  This 
improved  modeling  procedure  is  incorporated  into  an  inproved 
computational  algorithm  for  the  entire  flew  field,  yielding  an 
improved  modeling  and  computational  technique  that  is  unique  to  the 
present  research.  The  modeling  and  computational  technique  is 
validated  by  conparison  of  computed  results  with  available 
experimental  data.  Based  on  the  comparisons  a  favorable  assessment 
can  be  made  as  to  the  feasability  of  the  present  ccnputational 
procedure  for  routine  engineering  calculations.  Also,  it  is 
expected  that  the  calculated  heat  transfer  coefficients  presented 
may  in  themselves  be  useful  to  the  designer. 
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1.3  Brief  Description  of  Contents 

Hie  discussion  of  the  present  research  proceeds  as  follows. 
In  Chapter  II  the  state  of  the  art  of  investigation  of  turbulent 
recirculating  flow  is  reviewed,  with  enphasis  on  heat  transfer. 
Relevant  experimental  work  is  mentioned  and  attention  is  brought  to 
sane  of  the  most  common  closure  techniques  for  modeling.  Appendix 
A  discusses  the  mathematical  formulations  of  sore  of  these  closure 
techniques.  In  Chapter  III  the  mathematical  equations  actually 
used  for  the  present  computations  are  presented  with  the  detailed 
derivation  of  the  k-e  viscosity  model,  and  a  discussion  of  the 
present.  Ajuproveuents,  detailed  derivation  of  the  improvsnents  given 
in  Appendix  (B).  Chapter  III  also  includes  a  discussion  of  the 
less  successful  modeling  procedures  attempted  in  the  initial  stages 
of  this  study.  In  Chapter  IV  the  computat  ional  procedures 
including  the  treatment  of  the  near-wall  conditions  are  discussed. 
In  Chapter  v,  the  results  of  calculations  for  both  the  back-step 
and  pipe  expansion,  geometries  are  presented.  Comparison  is  made 
with  available  experimental  data  and  another  numerical  calculation 
(for  the  pipe  expansion) .  Heat  transfer  coefficient  correlations 
are  also  presented.  Finally,  in  Chapter  VI  conclusions  are 
discussed  concerning  the  calculation  procedures  and  their  general 
applicability,  and  finally,  suggestions  for  further  studies  are 
presented. 


CHAPTER  II 


REVIEW  OF  LITERATURE 

2.1  Review  of  Recirculating  Flew  Measurements 

Reliable  experimental  data  are  essential  for  the  development 
and  testing  of  turbulence  models.  This  is  simply,  because  the 
models  are  not  wholly  analytic;  for  instance,  they  contain  seme 
unknown  empirical  coefficients  (constants)  which  can  only  be 
estimated  by  comparison  to  experimental  data.  Hence,  availability 
of  accurate  measurements  means  that  the  turbulence  models  can  be 
adjusted  to  yield  reliable  predictions.  In  this  survey  of 
recirculating  flow  measurements  it  is  convenient  to  divide  the 
review  into  two  categories:  hydrodynamic  and  scalar  (temperature) 
flow  field  measurements. 

2.1.1  Hydrodynamic  Flow  Field  Measurements 

Measurements  of  hydrodynamic  flew  fields  for 
recirculating  shear  flews  are  much  more  extensive  them  those  of 
scalar  flow  fields.  Until  recently  measurement  of  turbulent  flow 
quantities  (even  the  time-averaged  quantities)  for  simple 
recirculating  flows,  such  as  the  back-step  was  very  difficult  [20] . 
The  flow  is  highly  turbulent  and  the  flow  reversals  make  hot  -  wire 
anemometers  generally  unsuitable.  Since  the  implementation  of 
laser-Dcppler  and  pulsed-vire  ananoroe  tries,  the  quality  of 
available  data  has  begun  to  iitprove.  The  following  literature 
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survey  represents  the  major  reported  research  on  reattaching  flows, 
with  emphasis  on  the  back-step  geometry.  A  similar  survey  was 
recently  compiled  by  Eaton  and  Johnston  120] . 

The  measurements  of  Hsu  (1950)  [21]  are  apparently  the  first 
experiments  to  study  both  mean  flow  and  turbulent  profiles  for  low 
speed  flows  over  a  large  step.  However,  his  turbulence  data, 
obtained  from  hot-wire  anemometers,  differs  significantly  from  many 
of  the  more  recent  experiments.  Another  early  study  is  that  of 
Abbott  and  Kline  [10]  in  1961.  They  studied  both  single  and  double 
backward-facing  steps  of  various  heights.  By  using  water  as  the 
working  fluid  and  employing  flow  visualization  techniques,  they 
obtained  detailed  measurements  of  mean  velocity  profiles  and 
reattachment  lengths.  This  work  includes  only  a  few  turbulent 
quantity  measurements  by  hot-film  anemcmetry.  Also  in  !961  Tani 
[22]  used  hot-wires  to  measure  mean  velocity,  turbulence  intensity 
and  turbulent  shear  stress  for  various  step  heights  for  both 
laminar  and  turbulent  incident  boundary  layers.  These  measurements 
apparently  represent  the  first  reliable  turbulence  data  for  the 
back-step. 

In  sane  of  the  more  recent  literature,  Bradshaw  and  Wong  [23] 
obtained  turbulence  measurements  in  the  entire  flow  region 
downstream  of  the  step.  These  data  were  then  used  to  obtain  sane 
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indirect  information  about  the  turbulence  structure  by  "numerical 
experiments" ,  i.e.,  by  forcing  empirical  calculations  to  match  the 
experiments.  This  then  enabled  them  to  study  in  more  detail  the 
shear-layer  mass  flew  recirculating  zone,  they  also  studied  the 
redeveloping  boundary  layer  and  noticed  that  for  at  least  52  step 
heights  downstream  the  attached  boundary  layer  characteristics 
("law  of  the  wall"  behavior)  were  not  fully  recovered. 

In  1975  Bradshaw  [24]  reviewed  and  identified  sane  complex 
turbulent  flews.  One  of  the  complex  flows  identified  was  separated 
flew.  Tto  contribute  to  the  understanding  of  the  separation  and 
redevelopment  phenomena,  recently  Chandrsuia  and  Bradshaw  [25] 
continued  the  experiments  of  Bradshaw  and  Wong  123] .  Their 
back-step  experiments  utilized  a  thin  laminar  approach  boundary 
layer,  mainly  bacause  in  the  reattached  region  full  development  is 
reached  much  faster  for  a  laminar  approach  boundary  layer  than  for 
a  turbulent  one  [26] .  Measurements  reported  include  mean  velocity 
and  Reynolds  stresses  (including  higher  order  correlations).  They 
concluded  because  of  the  inportance  of  the  triple  products,  any 
successful  modeling  procedure  should  reflect  these  effects  together 
with  the  'wall-effects'  in  the  pressure  -  strain  term.  The 
modeling  procedure  in  the  present  research,  to  be  discussed  later, 
does  include  seme  of  these  effects. 
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Narayanan  et  al.  [27]  have  measured  the  wall  static  pressure 
behind  steps  (of  varying  heights)  in  a  wind  tunnel.  Their  results 
indicate  reasonably  well  defined  similarity  patterns  in  the 
pressure  distributions  of  the  separated  flew  regions,  from  these 
distributions  they  were  able  to  estimate  reattachment  lengths. 

Kim  et  al.  [28]  measured  mean  velocity,  Reynolds  stresses  and 
intermittency  (fraction  of  time  for  which  the  flow  at  a  given  point 
is  turbulent)  in  the  entire  flow  field  of  the  back-step.  They 
noticed  that  downstream  of  rattachment  the  flaw  very  slowly  returns 
to  the  structure  of  the  attached  boundary  layer. 

For  the  double  facing  back-step  geometry,  Mehta  [29]  presented 
hot-wire  anemometer  results  of  mean  pressure  and  velocity  fields, 
turbulence  intensities  and  Reynolds  shear  stresses. 

Recently,  Moss  and  Baker  19]  presented  very  detailed 
measurements  by  using  the  plused-wire  anemometer  for  flews  over  a 
rearward-facing  step,  over  a  front-facing  step  and  over  a 
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profiles  and  turbulence  intensities. 

Laser  anemaretry  has  also  been  employed  in  the  study  of 
recirculating  flews.  Denham  and  Patrick  [30]  verified  the  use  of 
this  technique  by  presenting  detailed  velocity  profiles  for  laminar 
flows  beyond  a  back-step  in  a  plane  duct.  Etheridge  and  Kemp  [31] 
used  a  frequency-shifted  laser  to  study  the  flow  behind  the 
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back-step  in  a  water  channel.  They  presented  mean  velocities, 
turbulent  intensities  and  Reynolds  stresses.  By  also  using  the 
laser,  Smyth  L32]  studied  flws  over  a  double  backward-facing  step. 
Streamwise,  transverse  and  cross-stream  carnponents  of  velocity 
fluctuations,  turbulent  kinetic  energy  and  Reynolds  shear  stresses 
were  presented.  Armaly  et  al.  133]  also  presented  laser — Doppler 
measurements  of  velocity  distributions  and  reattactment  lengths  for 
the  back-step.  Results  were  for  laminar,  t  reins  it  ional  and 
turbulent  flow  of  air.  Bor  a  related  geometry,  Fraser  and  Siddig 
[11]  also  applied  laser  anerranetry  to  study  recirculation  due  to 
air  flow  over  a  normal  plane  wall  in  a  duct.  They  presented  mean 
and  fluctuating  velocities. 

For  the  pipe  expansion  Back  and  Roschke  [34]  measured  velocity 
profiles  for  laminar,  transitional  and  turbulent  approach  boundary 
layers,  the  fluid  being  water. 

All  of  the  above  studies  reported  reattachment  lengths  varying 
from  about  5.5  to  8.0  st  _  heights  downstream  of  separation. 
Reattachment  length  is  a  very  significant  quantity  in  recirculating 
flews  and  it  is  vorthw,'ile  to  give  sane  attention  to  the  factors 
that  affect  it.  For  example ,  changes  in  pressure  gradient  can 
cause  large  changes  in  the  reattachnent  length.  Recently,  Kuehn 
[35]  studied  this  effect.  Be  measured  the  influence  of  adverse 
pressure  gradient  on  resttachment  length.  The  resi'i  is 


are  shewn  in 
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Figure  (2),  which  indicates  an  increase  in  streanwise  pressure 
gradient  causes  an  increase  in  reattachment  length..  Narayanan  et 
al.  [27]  examined  the  effects  of  approach  momentum  boundary  layer 
thickness  on  reattachment  length.  Figure  (3)  summarizes  their 
experimental  results.  As  can  be  seen  the  effects  are  slight, 
particularly  in  the  range  5u/h  <1.5  (<5U  is  the  approach  boundary 
layer  thickness  and  H  the  step  height).  Finally,  the  effects  of 
approach  free  -  stream  turbulence  have  been  examined  by  Etheridge 
and  Kemp  [31] .  They  noticed  that  increasing  the  level  of  free 
stream  turbulence  results  in  a  decrease  in  reattachment. 

2.1.2  Scalar  Flow  Field  Measurements 

Mean  scalar  field  measurements  are  extremely  important  in 
understanding  the  fluid  mechanics  and  are  directly  applicable  to 
practical  situations.  Yet,  as  mentioned  earlier,  turbulence 
measurements  of  scalar  fields  are  very  rare  in  the  literature. 
Often  only  the  mean  scalar  field  is  reported;  usually  no  study  is 
n  ide  of  the  heat  transfer  coefficients,  scalar  fluxes  or  scalar 
fluctuations.  Conversely,  quite  a  few  heat  transfer  measurements 
are  concerned  only  with  scalar  measurments;  no  hydrodynamic  field 
measurements  are  included. 

Convective  recirculating  heat  transfer  flows  have  been 
reviewed  earlier  by  Hansen  and  Richardson  [36] ,  Fletcher  et  al. 
[37]  and  by  Aung  and  Watkins  [38] .  Aung  and  Goldstein  [6]  reported 


what  is  apparently  the  only  existing  low  speed  turbulent  heat 
transfer  data  for  the  back-step.  The  measurements  were  obtained  by 
using  a  Mach-Zehnder  interferometer.  Tenperature  profiles  and  heat 
transfer  coefficients  are  presented  for  the  entire  flew-  field. 
Much  earlier  than  Aung  and  Goldstein's  [8]  study,  Seban  et  al.  [6] 
studied  the  heat  transfer  over  the  back  step  for  high  speed  flews, 
with  the  upstream  plate  unheated.  For  the  pipe  expansion  Zemanick 
and  Dougall  139]  reported  what  are  also  the  most  complete  heat 
transfer  measurements  for  the  entire  flow  field.  Also  in  related 
work,  much  earlier  than  this  study,  were  the  measurements  of  Ede  et 
al.  [40] «  Tn  their  study  they  did  not  enphasize  heat  transfer  in 
the  redeveloping  region.  Recently,  Sparrow  and  O'Brien  [41] 
presented  heat  transfer  data  for  the  heated  side  wall  in  a  pipe 
expansion  geometry.  The  heat  transfer  coefficients  were  determined 
by  a  napthalene  sublimation  technique.  It  should  be  reemphasised 
that  it  is  unfortunate  than  none  of  the  above  studies  measured  the 
hydrodynamic  flow  field  in  conjunction  with  the  heat  transfer. 

Some  heat  transfer  measurements  for  related  geometries  are:  a 
double  step  by  Seki  et  al.  [42,43],  a  blockage  in  a  pipe  by  Koram 
and  Sparrow  [44] ,  flow  over  a  disc  by  Srayth  [45] ,  high  speed  flow 
over  a  cavity  by  Haugen  and  Dhanak  [46]  and  flow  over  blunt  plate 
by  Ota  and  Kon  [47,48]. 

All  of  the  recirculating  flow  heat  transfer  measurements 
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confirm  a  peak  in  the  surface  heat  transfer  coefficient  near  the 
reattachment  point,  the  peak  value  being  mary  times  larger  than  the 
fully  developed  (or  attached)  value.  Correlations  of  experimental 
heat  transfer  data  are  very  rare  in  the  literature,  probably 
because  of  the  wide  range  of  discrepancies  between  the 
measurements.  However,  for  the  back-step  arva  for  supersonic  flews, 
a  correlation  procedure  was  attenpted  by  Lamb  [49] . 
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2.2  Turbulent  Modeling  For  Recirculating  Flows 

A  great  deal  of  modeling  success  has  been  achieved  for 
turbulent  boundary  layer  flew,  with  and  without  heat  transfer.  For 
these  flows  empirical  relations  such  as  the  'law  of  the  wall'  were 
adequate.  Such  approaches  fail  for  turbulent  separated  flows, 
mis  then  led  to  the  development  of  other  theoretical  models,  both 
simple  and  complex.  The  advantages  of  the  simple  models  are  that 
they  enable  the  development  of  very  versatile  and  useful  numerical 
codes.  Employing  these  cedes  furnished  and  emphasized  the  need  for 
more  realistic  and  advanced  turbulent  models. 

As  seen  in  Chapter  I,  in  order  to  solve  the  mean  momentum  and 
mean  scalar  equations,  Bgns.  (1.2.3)  and  (1.2.5),  one  needs 
information  on  the  correlations,  u^u!  and  6* .  Before  reviewing 
the  modeling  procedures  employed  by  researchers  for  these  terms  in 
heat  transfer  recirculating  flaws,  a  generalized  discussion  of  seme 
of  the  proposals  for  turbulent  closures  is  presented. 

2.2.1  Sane  Closure  Techniques 

mere  are  a  few  extensive  reviews  in  the  literature  that 
examine  the  various  closure  techniques.  For  example,  four  such 
reviews  are:  Bradshaw  [50] ,  Rodi  [51] ,  Reynolds  [52]  and  Mellor  and 
Herring  [53] .  Modeling  is  better  accaiplished  with  an  insight  into 
the  physics  of  turbulence.  Such  an  insight  is  given  by  Bradshaw 
[50] .  In  this  section  modeling  of  both  hydrodynamic  and  scalar 
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flew  fields  is  discussed.  Appendix  A  contains  the  mathematical 
expressions  of  the  models  discussed  below. 

(a)  Eddy  Viscosity  Models 

This  closure  procedure  is  based  on  Boussinesq’s 
(1877)  [54]  idea  of  a  "turbulent  visocity" .  He  suggested  that  the 
local  turbulent  shear  stress  be  related  to  the  eddy  viscosity  and 
the  local  mean  velocity  gradient,  i.e.,  analagous  to  laminar 
stresses.  The  turbulent  (eddy)  viscosity  is  assumed  a  property  of 
the  flew  field  and  not  a  fluid  property.  For  scalar  fluxes  use  is 
made  of  the  ratio  of  eddy  viscosity  to  eddy  diffusivity,  denoted  by 
the  turbulent  Prandtl  number.  Closure  is  accomplished  by 
specifying  the  eddy  viscosity. 

Two  common  procedures  are  in  use,  one  expresses  the 
eddy  viscosity  algebraically  in  terms  of  mean  quantities  and  the 
other  fcy  differential  equations.  The  simplest  algebraic  expression 
is  the  mixing  length  hypothesis  of  Prandtl  [55J .  An  improvement  of 
this  length  scale  was  suggested  by  Vein  Driest  [56] . 

The  first  level  of  improvement  of  the  algebraic 
viscosity  model  is  the  so  called  one-equation  model.  In  this 
model,  in  approximating  the  turbulent  viscosity,  a  differential 
equation  for  a  velocity  scale  (in  terms  of  the  turbulent  kinetic 
energy,  k)  is  introduced  together  with  the  mixing  length  scale. 
For  boundary  layer  flews  this  procedure  was  tested  by  Bradshaw  et 
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al. [57] .  The  turbulent  kinetic  energy  is  an  important  parameter, 
for  it  distinguishes  turbulent  frcm  laminar  flew  and  is  well 
documented  by  experiments. 

A  major  problem  in  turbulent  modeling  techniques  is 
the  choice  of  necessary  and  adequate  equations.  Most  work  features 
the  diffusion  of  kinetic  energy,  e,  as  a  key  quantity.  An  improved 
modeling  procedure  involves  approximating  the  viscosity  in  terms  of 
differential  equations  for  the  velocity  and  length  scales  (done  in 
terms  of  k  and  e  ) .  This  procedure  is  commonly  referred  to  as  the 
two  -  equation  model  of  turbulence  or  simply  the  k-e  model.  The 
standard  k-e  model  (high  Reynolds  number)  requires  specification  of 
'wall-functions'  [58],  Same  improvements  of  the  standard 
•wall-functions’  are  given  by  Chieng  and  Launder  [59] .  Jones  and 
Launder  160]  extended  the  standard  k-  e  model  to  include  near-wall 
effects;  no  specification  of  the  'wall-functions'  outside  the 
viscous  sublayer  is  needed.  Because  of  the  simplicity  and  success 
of  the  k-  e  model  many  recent  studies  were  focused  on  inproving  the 
deficiencies  of  the  model.  Ihe  formulation  of  the  e  differential 
equation  has  been  the  topic  of  much  controversy.  Re-examining  the 
equation,  Hanjalic  and  Launder  [61,62]  'sensitized'  the  diffusion 
of  energy  effects  by  considering  the  importance  of  near-wall 
turbulence  and  the  influence  of  normal  strains.  Another  drawback 
of  the  model  is  that  it  uses  an  enpirical  constant  in  the 
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formulation  of  the  turbulent  viscosity.  Recently,  Ljuboja  and  Rodi 
[63] ,  Leschziner  and  Rodi  [64] ,  Humphrey  and  Pourahmadi  [65]  and 
Gooray,  Watkins  and  Aung  [66]  replaced  these  constants  by  a 
function  derived  frcm  an  algebraic  stress  models  The  improvement 
procedure  is  outlined  in  a  later  chapter.  In  the  present  study  an 
analagous  calculation  was  performed  for  the  turbulent  Prandtl 
number  and  is  discussed  in  a  later  section.  An  earlier  study  by  So 
[67] ,  attempted  a  similar  procedure,  but  by  using  the  curvature 
modifications  of  Bradshaw  [68] . 

(b)  Stress  and  Flux  Transport  Models 

Rather  than  assume  a  direct  relation  between  the 
shear  stress  and  mean  strain  field  (Boussinesq's  idea),  transport 
equations  for  the  Reynolds  stresses  and  scalar  fluxes  are 
developed.  The  advantages  of  this  procedure  are  that  solutions 
start  at  the  wall  itself,  thus  no  specification  of 
'wall-functions',  or  assumptions  as  to  the  turbulent  Prandtl  number 
are  required,  and  it  automatically  includes  such  effects  as 
streamline  curvatures  and  pressure-strain  (scalar)  interactions . 
However,  it  does  require  solving  some  5  to  7  transport  equations  of 
turbulence. 

This  modeling  procedure  is  much  more  reliable  for 
predicting  non-equilibruim  flow  (which  include  recirculating  and 
mixing  flows),  i.e.,  when  using  the  turbulent  viscosity  procedure. 
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the  forcing  of  the  shear  stress  to  instantaneously  respond  to 
changes  in  the  mean  strain  field  is  not  conceptually  valid  in  such 
flows. 

•Thus,  the  closure  of  the  approximate  differential 
equations  for  the  Reynolds  stresses  (u  uju)  and  scalar  fluxes  (uie ' ) 
offer’s  a  very  premising  class  of  turbulence  models  for  the  solution 
of  complicated  turbulent  flews.  Hie  development  of  the  closure 
procedures,  including  procedures  for  scalar  fluxes,  can  ba  found  in 
many  references,  e.g.  Launder,  Reeoe  and  Rodi  [17] ,  Launder  [69]  , 
Hanjalic  and  Launder  [71]  (extension  of  the  model  to  include 
near-wall  effects) .  Closure  of  the  stresses  and  fluxes  requires 
approximations  for  the  pressure  -  strain/scalar  and  triple 
correlation  terras.  The  various  approximations  suggested  for  these 
terns  differ  significantly;  for  this  reason  much  present  research 
is  still  directed  in  this  area.  For  example,  Gormack  et  al.  [72] 
developed  an  improvement  of  the  Hanjalic  and  Launder  [71]  triple 
velocity  correlation  closure.  However,  because  of  the  model’s 
complexity  Gormack  et  al.  [72]  reoenntanded  the  one-parameter  model 
of  Hanjalic  and  Launder  [71]  for  meet  basic  flew  situations.  This 
conclusion  was  arrived  at  after  .making  comparisons  with  various 
other  models  and  experiments  for  axl symmetric  channel  flow,  pipe 
flew,  wall-jets  and  mixing  layers.  Recently,  both  So  [73]  and 
Spesiale  [74]  indepervlently  further  studied  the  pressure-strain 
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modeling  of  Hanjalic  and  Launder  [71] .  In  his  work,  So  [73] 
included  terms  which  account  for  rapid  distortion  of  the  shear 
stress  and  the  return  to  isotropy  process.  However,  the 
improvements  were  found  to  be  minor  when  compared  to  the  Rotta 
return  to  isotropy  model  used  by  Hanjalic  and  Launder  [71] .  This 
indicated  that  much  improvement  of  the  pressure-strain  terra  is  not 
possible  by  just  refining  the  simple  Stotta  model,  a  procedure 
suggested  by  Launder  [69] .  With  this  in  mind  Speziale  [74] 
redefined  the  pressure  ■  strain  terra.  Proving  this  terra  to  be 
invariant,  he  developed  a  representation  for  it  that  is  dependent 
on  the  Reynolds  stresses,  mean  velocity  gradients  and  a  master 
length  scale.  However,  this  inprovesnent  still  has  to  be  tested  to 
determine  its  usefulness  and  range  of  validity.  The  above  modeling 
procedures  do  not  correctly  predict  the  complex  turbulent  scales 
that  are  present  in  corner  flews.  Recently,  Gessner  and  Emery  and 
Gessner  [75,76]  improved  the  Launder,  Recce  and  Rodi's  model  [71] 
to  better  simulate  corner  flows. 

Modeling  of  the  scalar  fuxes  have  been  carried  out 
by  Launder  [69].  Gibson  and  Launder  [77]  extended  the  models  to 
include  buoyancy  effects  for  thin  shear  flews  remote  frau  wed. Is. 
Recently,  Samara weera  [78]  performed  extensive  testing  of  the 
various  scalar  flux  models  for  boundary-layer  type  flews. 

Analogous  to  the  pressure-strain  modeling  is  the 


30 


pressure-  scalar  modeling.  This  procedure  is  outlined  in  Gibson 
[79] .  The  empirical  constants  appearing  in  the  pressure-scalar 
model  can  be  obtained  frcm  the  data  of  Antonia  and  Danh  [80] . 

In  an  effort  to  reduce  the  number  of  transport 
equations  to  be  solved,  sirrpler  versions  of  the  Reynolds  stresses 
and  scalar  fluxes  closure  schemes  have  been  derived.  One  such 
simplification  is  the  so  called  algebraic  stress  closure.  Hie 
mathematics  of  the  procedure  is  given  in  the  next  chapter,  and  is 
similar  to  that  formulated  by  Launder  [69] .  Gessner  and  Emery  [77] 
and  Gessner  and  Po  [81]  applied  this  model  to  flow  in  a  rectangular 
duct,  and  found  that  the  model  is  very  accurate  when  only  two 
empirical  constants  and  a  generalized  mixing  length  are  specified. 
Also,  So  [82]  showed  that  this  model  can  be  used  to  predict  the 
behavior  of  different  kinds  of  shear  flews,  i.e.,  flows  over  curved 
surfaces.  Algebraic  modeling  of  the  scalar  flux  equation  has  been 
done  by  Gibson  [79]  and  Gibson  and  Launder  [77] . 

As  previously  mentioned,  algebraic  modeling  has  been 
used  to  estimate  the  constant  in  the  eddy  viscosity  model,  thereby 
simultaneously  accounting  for  streamline  curvature  and  pressure  - 
strain  effects.  Bradshaw  [68]  reported  an  early  study  on  the 
importance  of  these  effects.  Irwin  and  Smith  [83]  used  the 
algebraic  stress  model  to  also  demonstrate  these  effects.  However, 
because  of  the  difficulties  involved  in  specifying  length  scales 
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and  an  eirpirical  constant  in  Bradshaw's  168]  model,  the  algebraic 
stress  nodeling  seems  more  adaptable  to  improving  the  k-e  model. 
2.2.2  Application  to  Recirculating  Flows 

The  mixing  length  models  can  be  used  for  satisfactory 
prediction  of  two-dimensional  boundary  layer  flows.  Since  these 
flows  are  self  preserving  and  have  local  equilibria im  regions,  the 
single  length  scale,  with  correctly  chosen  empirical  constants  can 
yield  satisfactory  results.  However,  applications  to  recirculating 
flows  have  only  met  moderate  success  at  best  [91] . 

Use  of  the  two-equation  (k-f.)  model  is  an  obvious 

inprovement,  since  it  incorporates  the  transport  effects  of  the 
length  scale  into  the  turbulence  model.  A  review  of  the 

applicability  and  validity  of  this  model  is  demonstrated  by  Launder 
and  Spalding  [58]  for  many  flow  situations.  Included  is  flow 
recirculation,  i.e.,  a  film  cooling  application  and  flow  in  a 
cavity. 

After  the  inception  of  the  present  study,  a  few 

researchers  have  attenpted  modeling  recirculating  flows  for  the 

back-step  and  pipe  expansion.  In  all  cases  known  to  the  author, 
the  k-r  model  was  used.  For  momentum  transfer  studies,  Taylor  et 
al.  [84]  incorporated  this  model  with  a  finite  element  solution 
procedure  to  yield  velocity  profiles,  turbulent  intensity  and 
reattachment  lengths.  Oliver  [85]  studied  heat  transfer  over  the 
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back-step  and  compared  his  results  with  the  experiments  of  Seban 
[6]  for  high  speed  flows.  His  predictions  overestimate  the  heat 
transfer  data  in  the  entire  flow  field,  particularly  in  the 
reattachment  region.  Chieng  and  Launder  [59]  studied  the  pipe 
expansion  geometry',  by  using  both  the  standard  k-e  model  and  the 
low  Reynolds  number  version  of  the  model.  They  caipared  heat 
transfer  coefficients  with  the  data  of  Zenamick  and  Dougall  [39] . 
A  more  detailed  discussion  of  the  Chieng  and  Launder  [59] 
calculations  is  presented  later. 

An  early  analytic  analysis  of  heat  transfer  for  turbulent 
separated  flows  was  performed  by  Spalding  [S6j -  He  considered  a 
one  dimensional  analysis  near  the  wall  and  postulated  relationships 
for  the  generation,  diffusion  and  dissipation  of  turbulent  energy 
in  terms  of  the  turbulent  kinetic  energy  and  empirical  constants. 
Solving  the  differential  equation  for  k,  the  Stanton  number 
dependence  was  obtained.  However,  as  Spalding  [86]  suggested  a 
two-dimensional  analysis  with  a  more  realistic  turbulent  model  is 
needed  so  as  to  quantitatively  determine  velocity  and  temperature 
profiles  and  turbulent  intensities. 

Modeling  procedures  for  recirculating  free  shear  flows 
(such  as  jets)  have  also  been  studied.  Such  flows  have  been 
reviewed  by  Durst  and  Rastogi  [87]  and  by  Leschziner  and  Rodi  [64] . 
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2.3.  Review  of  Computational  Methods  for  Recirculating  Flows 

Conputation  of  recirculating  flews  has  been  a  major  problem 
for  decades  new.  For  laminar  flows  the  hydrodynamic  flow  field  for 
the  back-step  was  reliably  predicted  by  the  vorticity/stream 
function  formulation  of  Roache  and  Mueller  188] .  Turbulent  flow 
field  simulations  of  recirculating  flow  by  solving  the  complete 
Navier-StcJces  equations  have  increased  with  the  advent  of  high 
speed  computers.  As  discussed  in  the  previous  section  great 
difficulties  exist  in  the  formulation  of  a  set  of  differential 
equations  to  model  the  flow;  and  as  expected,  numerical 
canpu  tat  ions  of  these  equations  are  also  very  tedious.  For  this 
reason  many  investigators  have  concentrated  their  efforts  at 
locking  at  the  numerical  schemes.  A  review  of  seme  of  this 
research  is  discussed  below. 

Recently,  Atkins  et  al.  189]  examined  various  numerical 
schemes  for  separated  flows  and  compared  the  different  results 
obtained  for  the  back-step  geometry.  They  concluded  that  for  low 
Reynolds  nunber  by  using  pure '  upwind  differencing  to  approximate 
the  convective  terms  yields  reattachment  lengths  about  8  percent 
lower  tiian  that  obtained  by  a  corresponding  central  differencing 
approximation  procedure.  At  high  Reynolds  number,  the  upwind 
difference  prediction  agrees  well  with  experimental  results,  but 
the  numerical  specification  of  the  inlet  condition  is  very 
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critical.  Also,  Minh  and  Chassaing  [90]  recently  studied  turbulent 
numerical  calculations;  for  both  boundary  layer  and  recirculating 
type  flows.  Emphasis  was  not  placed  on  turbulent,  modeling,  but  on 
choice  of  unknown  variables  ( vert icity /stream  function,  w  or 
pressure/ velocity,  u,  v,  p]  and  the  influence  of  inlet  conditions. 
They  found  that  the  (u  v  p)  primitive  variable  formulation  was 
better  than  the  {  )  method  for  yielding  pressure  fields  and 
turbulent  quantities  (k,  u'v').  Castro  [91]  also  noted  that  the 
accuracy  of  the  cannon  differencing  schemes  is  often  suspect.  He 
showed  that  upwind  differencing  is  not  accurate  near  the  sharp 
edges  at  which  the  flow  separates,  and  found  that  skew  differencing 
schemes  seem  to  be  better.  However,  skew  differencing  techniques 
have  not  been  fully  tested  and  verified  for  recirculating  flows. 
Richards  and  Crane  [92]  developed  a  new  central  differencing  scheme 
and  found  it  to  be  more  accurate  than  the  conventional  upwind 
differencing  scheme.  However,  the  scheme  will  only  yield 
convergence  if  there  is  an  upper  lhi.it  on  a  quantity  defined  as  the 
grid  Reynolds  number.  Unfortunately,  to  satisfy'  this  condition 
requires  an  excessively  fine  mesh  system  and  restricts  the  Reynolds 
number  to  extremely  low  values.  Biringen  [93]  solved  the  partial 
differential  equations  by  a  time  explicit  finite  differencing 
scheme,  a  modified  Lax-Wendroff  method.  He  found  that  values  of 
the  turbulent  kinetic  energy  were  very  sensitive  to  the  aipirical 
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input;  whereas,  mean  velocities  and  stresses  were  not  so  sensitive. 

A  recent  review  by  Barrett  and  Demunshi  [94]  examined  finite 
element,  finite  difference  and  finite  integral  methods  for 
recirculating  flews.  They  concluded  that  finite  difference 
procedures  seem  to  more  reliably  predict  the  recirculating  flow 
field. 

Initially  many  calculation  procedures  utilized  the  hybrid 
centred/upwind  scheme  for  the  convective  terns  and  centred 
differencing  for  the  diffusive  terms.  This  procedure  was  developed 
by  Patankar  [95] .  Very  recently,  Patankar  [96]  reviewed  the 
procedure.  The  procedure  suffers  fruu  inaccuracies  due  to 
truncation  error  and  false  diffusion.  In  light  of  this,  many 
recent  studies  have  focused  on  inproving  the  hybrid  scheme.  Two 
improvenent  procedures  were  studied  recently,  namely,  the  hybrid 
central/ skew-upwind  differencing  scheme  of  Raithby  [97]  and  the 
quadratic  differencing  scheme  of  Leonard  [98] .  The  quadratic 
scheme  is  incorporated  in  this  research  and  will  be  discussed  in  a 
later  chapter.  Practical  evaluation  of  these  schemes,  very 
recently.-  were  reported  by  Leschziner  [99] ,  Leschziner  and  Rod! 
[64]  an.  Han,  Hurtphrey  and  launder  [100] .  The  numerical  procedure 
enplcyed  in  this  study  employs  the  improvement  of  Leonard  [98]  on 
the  Patankar  procedure  [95] ,  as  developed  by  Gosman  [101] . 


2.4  SUMMARY 
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Presently,  the  information  on  the  complex  turbulent  flows  of 
the  backward  step  and  the  sudden  pipe  expansion  is  sufficient 
enough  for  qualitative  discussion  of  these  flews.  However, 
engineering  applications  require  quantitative  predictions  of  sane 
of  these  complex  flows.  In  the  foregoing  sections  such  prediction 
possibilities  were  explored  by  locking  at  the  existing  experimental 
measurements  and  modeling  procedures  for  turbulent  recirculating 
heat  transfer,  with  emphasis  on  the  back-step  and  pipe  expansion 
geometries.  The  existing  numerical  procedures  were  also  examined 
with  the  conclusion  being  that  they  are  wall  tested  and  quite 
adequate  for  use  in  studying  recirculating  flows.  So  with 
confidence  in  the  numerical  procedure,  the  validity  and 
applicability  of  the  turbulence  models  can  be  tested  to  same 
degree.  However,  for  such  modeling  to  be  more  meaningful  there  is 
a  need  for  more  extensive  heat  transfer  measurements  for  these 
basic  engineering  geore  tries. 

As  a  final  remark,  the  turbulent  modeling  procedures  do  not 
take  into  account  the  intermittent  nature  of  turbulence  that  simple 
experiments  can  detail.  This  is  due  to  the  fact  that  in  seeking  a 
nodel  for  tractable  ccqputer  solution,  resorting  to  the  averaging 
of  the  equations  of  motion  smears  out  this  time  dependency.  Such 
effects,  however,  are  probably  more  dominating  in  free  shear  flews 
than  in  wall  bounded  flows  and  moreover,  the  time  scales  involved 
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are  probably  riot  significant  when  considering  application  to  heat 
exchange  equipment.. 
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CHAPTER  III 


MAIHEIMTCAL  AND  PHYSICAL  MODEL  fcTJR  GEOMETRIES  CONSIDERED 
3rl  General  Remarks 

As  discussed  previously  in  Chapter  X,  solution  of  the 
Navier-Stokes  and  energy  equations  for  turbulent  flews,  Eqns. 
(1.2.3)  and  (1.2.5),  requires  closure  of  the  Reynolds  stresses  and 
scalar  fluxes.  In  Chapter  II  the  various  closure  schemes  were 
described.  In  the  present  study,  variations  of  one  such  closure 
scheme,  turbulent  viscosity  modeling  (modeling  based  on  the 
Boussinesq  [54]  approximation) ,  were  applied  to  two  recirculating 
flow  geometries  where  consistent  experimental  data  exists:  air  flow 
at  moderate  Reynolds  numbers  over  two-dimensional  rearward-facing 
steps  and  sudden  pipe  expansions.  Successful  predictions  of  the 
turbulent  flow  over  the  two-dimensional  rearward-facing  steps  and 
axisynmetric  pipe  expansions  were  eventually  obtained  by  solving 
numerically  the  governing  partial  differential  equations  as 
formulated  for  the  k re  turbulent  viscosity  model  by  Launder  and 
co-workers,  cf  [19,58,61].  Although  the  present  study  focuses  on 
the  so  called  turbulent  viscosity  models,  details  of  other 
turbulence  closure  schemes  as  well  as  details  of  the  turbulent 
viscosity  models  are  presented  in  Appendix  (A)  for  convenient 
reference  purposes. 
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3.2  The  Boussinesq  Approximation 

As  mentioned  before,  the  -  closure  procedure  used  in  the  present 
study  incorporates  the  Boussinesq  [54]  turbulent  viscosity 
hypothesis,  which  relates  the  Reynolds  stresses  to  mean  rates  of 
strains  by  a  turbulent  viscosity  concept.  According  to  this 
formulation,  the  various  ccrponents  of  the  Reynolds  stresses  are. 


-p  U  V 


pk 

pk 


(3.2.1) 

(3.2.2) 

(3.2.3) 


As  before,  is  the  turbulent  viscosity.  Assuming  similar 

relationships  for  the  heat  fluxes. 
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Where  I‘t  («  Ut/Prt)  is  the  turbulent  diffusLvity  of  temperature. 
The  turbulent  viscosity  has  to  be  modeled  and  1’^  is  approximated  by 
using  the  turbulent  Prandtl  number,  Pr^,  i.e.,  the  ratio  of 
turbulent  viscosity  to  turbulent  diffusivity  (ht/rt)‘ 
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Prt  »  Ut/Tt 

3.3  Fluid  and  Thermal  Transport  Equations 

By  employing  the  Boussinesq  [54]  approximations  of  the 
previous  section,  the  time  independent  mean  icomentisn  and  energy 
equations,  Bgns.  (1.2.1),  (1.2.3)  and  (1.2.4),  for  an 

incompressible  fluid  (for  each  component)  be  carve: 
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(3.3.2) 


(3.3.3) 


where,  j  =  0  applies  to  the  back  step  and 


41 

j  =  1  applies  to  the  pipe  expansion 
In  the  above 

Pgff  =  y  +  Pt/  where  is  the  turbulent  viscosity 

reff  •  **  +  ^t.  ,  where  Prt  is  the  turbulent  Prandtl 

Pr  Pr. 

^  number 

In  the  foregoing  equations,  for  notational  simplicity  the  usual 
overbars  have  been  emitted  for  the  mean  fluid-dynamic  and 
temperature  variables.  As  shown  in  Figure  (4) ,  y  is  measured  from 
the  free  stream  of  the  back-step  or  the  pipe  centerline,  and  x  is 
the  distance  along  the  plate  or  pipe  axis  fran  the  vertical  face. 

Solution  of  the  above  equations  requires  solving  a  set  of 
coupled  non-linear  elliptic  partial  differential  equations,  with  a 
provision  for  determination  of  the  pressure  field.  Initial 
conditions  are  discussed  in  a  later  section  and  boundary 
conditions,  together  with  the  solution  algorithm  are  discussed  in 
the  next  chapter. 

Specification  of  yj.  and  Pr^,  so  as  to  solve  the  system  of 
equations,  constitutes  a  major  part  of  this  research.  In  the 
initial  phase  of  this  study,  the  simplest  turbulent  viscosity  model 
(the  mixing  length  model)  together  with  a  constant  value  for  the 
turbulent  Prandtl  number  was  used.  As  will  be  discussed  later, 
because  of  the  presence  of  complexities  in  the  turbulence 
associated  with  recirculating  heat  transfer  flews,  the  mixing 
length  model  was  replaced  with  a  somewhat  more  sophisticated 
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turbulent  viscosity  the  k-e  model  which  gave  better  results.  As 
the  study  progressed,  it  became  evident  that  in  order  to  further 
improve  overall  predictions,  modifications  of  the  k-e model  were 
necessary.  Numerical  experiments  were  then  performed  to  determine 
the  most  reliable  form  of  the  k-  e  model  that  could  be  applied  for 
engineering  predictions  of  recirculating  heat  transfer.  Hie  final 
calculations  in  this  study  use  the  modified  standard  k-e  model, 
together  with  a  variation  of  this  model  appropriate  for  wall  shear 
flews.  Functional  relationships  were  derived  for  two  coefficients 
appearing  in  the  model,  a  coefficient  that  relates  the  viscosity  to 
the  turbulent  length  and  velocity  scales,  and  the  turbulent  Prandtl 
number. 

The  remainder  of  this  chapter  is  devoted  primarily  to  the 
formulations  of  the  various  viscosity  models  tested  in  this  study. 
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3.4  Turbulent  Prandtl  Number  and  Viscosity  Models 
3.4.1  Turbulent  Prandtl  Number 

In  most  previous  numerical  calculations  Prt  has  been 
assuned  to  be  a  constant  value,  e.g.,  for  air  a  value  of  0.9.  This 
value  has  been  substantiated  by  experiments  for  regions  outside  the 
viscous  sublayer  [102].  Hoover,  in  the  viscous  sublayer  as  shown 
by  Blackwell  [102]  the  value  is  significantly  different.  Crawford 
and  Kays  [103]  developed  an  empirical  expression  which  simulates 
the  data  very  well. 


+  c  Pet  (Prt00) 
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where  Pr^  is  the  free  stream  value  of  0.9,  c  »  0.2  and 


Pet  -  (Ut/u)  Pr  , 


Pr  (molecular)  *  0.7  (for  air) 


In  the  above,  the  turbulent  viscosity,  yt,  was  approximated  by 
using  the  mixing  length  model. 

By  using  a  Taylor  series  expansion  procedure,  Antonia 
U 14]  developed  a  simpler  analytic  approximation  for  the  behavior 


of  Prt  in  the  Buffer  zone  (5  <  y  +  <  20). 

Pr*  s  0.61  £1-0. 045V4-) 

where  y4-  is  the  non-dimensional  normal  co-ordinate, 
y4"  »  y  u^/  v  ,  a,.-  (t^p)3*,^  is  the  wall  shear  stress. 


(3.4.2) 
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In  the  present  st  iy  preliminary  calculations  used  Bqn. 
(3.4.2)  in  the  viscous  sublayer  and  a  value  of  0.9  elsewhere.  Such 
a  procedure  is  somewhat  inconsistent  for  the  recirculating  region, 
since  Eqn.  (3.4.2)  was  derived  for  a  developed  boundary  layer.  In 
subsequent  calculations,  the  above  approach  was  abandoned  for  a 
more  generalized  description  of  Pr^  based  on  algebraic  stress 
modeling.  Since  the  ooncept  and  derivation  of  this  new  expression 
is  similar  to  that  of  the  Improved  turbulent  viscosity  model, 
discussion  of  it  is  deferred  to  Section  (3.4.4)  where  both  are 


discussed  concurrently. 


3.4.2  Mixing  Length  Model 

The  mixing  length  hypothesis  is  the  simplest  closure 
procedure.  In  this  model  the  viscosity  is  approximated  by  a  single 
length  scale.  During  the  initial  phase  of  the  present  study 
numerical  experiments  were  performed  with  a  modified  Van  Driest 
[105]  length  scale,  where, 


in  which 

*o  “  K  y  [l-exp  ( -y+V A+)J 

where  A+  is  normally  taken  to  be  a  constant  value  of  26, 


(3.4.3) 


(3.4.4) 
to  conform 


with  the  fully  turbulent  profile  remote  from  a  wall.  Initial 
calculations  were  performed  incorporating  the  above  length  scale, 
but  with  A+  given  a  functional  behavior,  as  discussed  by  Crawford 
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and  Kays  [103.]  Unfortunately,  the  heat  transfer  predictions  by 
this  simple  model  grossly  underestimated  the  data,  particularly  in 
the  recirculating  region.  Such  a  result  however,  is  anticipated, 
for,  in  the  recirculating  (ana  redeveloping )  region  it  should  not 
be  expected  that  a  single  length  scale  will  accurately  predict  the 
conplf’X  turbulence  phenomena  that  exist.  After  the  unsuccessful 
initial  numerical  experiments,  the  mixing  length  approach  was 
abandoned  in  favor  of  the  k-e  model. 

3.4.3  k-e  Model 

The  k-  e  viscosity  approximation  procedure  is  one  of  the 
most  successful  and  widely  used  closure  techniques.  The  standard 
approximation  for  the  viscosity  by  two  scales  is, 

1/2  k  3/2  ' 

yt  “  Cy  P  k  ~e —  (3.4.5) 

where  C  ^  is  normally  assured  a  constant,  k^  (k  is  the  turbulent 
kinetic  energy,  i.e.,  k  *  u^u!/2;  is  the  vel  ^  scale  and  k  3/2/ e 
is  the  length  scale,  e  being  the  dissipation  rate  of  turbulence 
energy.  Transport  differential  equations  for  k  and  e  are  required 
and  will  be  presented  shortly  with  non-essential  details  given  in 
Appendix  a.  Hcwever,  first  scras  remarks  on  are  appropriate. 

Early  calculations  in  the  present  research  1106]  used  the 
established  constant  value  of  0.09  for  [58] ,  but  subsequent 
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studies  indicated  the  need  to  improve  this  model  through  giving 
a  functional  value.  The  rationale  for  this  improvement  is  that  in 
recirculating  flews  the  effects  of  streamline  curvatures  and 
pressure-strain  interactions  are  significant  factors  that  influence 
the  recirculating  process  [63,64,651.  Neglecting  these  effects 
results  in  a  shorter  recirculating  zone.  In  the  k-e  model  the  most 
consistent  manner  in  which  both  these  effects  can  be  included  is 
through  .  As  in  the  case  of  Pr^  a  generalized  expression  based 
on  algebraic  stress  modeling  can  be  derived  and  further  discussion 
is  deferred  to  Section  (3.4.4). 

In  this  study  two  forms  of  the  k-e  model  are  used,  the 
standard  or  "  h  igh-Reynold  s-number"  version  158]  and  the 
lcw-Reynolds-number  form  proposed  by  Jones  and  launder  1106] .  In 
the  generalized  equations  presented  below,  terms  appearing  within 
the  broken  lines  are  emitted  when  using  the  standard  version  (Eqn. 
3.4.5)  of  the  model.  The  effective  viscosity  and  diffusivity  terms 
appearing  in  Eqns.  (3.3.2)  -  (3,3.4)  are  given  by  (approximations 
for  Prt  were  discussed  previously) , 

^ff"u  +  ^t  ; 

r  +  ^ 

eff  ^Fr  Prt 
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The  transport  of  turbulent  kinetic  energy,  k,  is  expressed  by 
(Appendix  A) 


similarly  be  expressed  by,- 

p  —  (ue)  +  p  (ve)  m  —  T  (  '  U  +' 

3x  3y  jr  >  3x  *•  '  i  1 


In  the  preceding  equations, 

G52[(!i)2+(f)2  +  H!r)2]  +  (l7  +  fe)2 

C^f  Cei,  C$2'  ak  ae  are  cons  tints  of  the  standard  model  and  are 
presented  in  Table  (3.1).  The  low-turbulent  Reynolds  number 
functions  are  taken  as, 
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where  values  for  and  Ar  are  also  presented  in  Table  (3.1). 

•ttie  standard  k-e  model  as  originally  described  by  Launder  and 
Spalding  [19]  does  not  have  the  capability  of  including  details  of 
flow  processes  involving  the  molecular  viscous  effects  near  a  wall. 
Away  from  a  wall  ( in  the  logarithmic  region) ,  where  the  influences 
of  molecular  viscosity  are  insignificant  (y  <  yt),  che  standard 
eqi,rhlCTis  for  k  and  e  are  valid.  Thus  to  account  for  the  viscous 
diffusion  of  k  and  e  close  to  the  wall.  Launder  and  Jones  [60] 
added  the  molecular  viscosity  dependent  terns  to  the  standard 
equations.  The  viscosity  dependent  functions  are  related  to  the 
turbulent  Reynolds  number  R^pk^/ey) ,  which  has  a  dominating 
influence  near  the  wall.  These  additional  terms  are  based  on 
empirical  relationships  obtained  fran  boundary  layer  type  flews;  as 
such,  they  are  not  based  on  exact  mathematical  formulations. 

Apart  from  the  lcw-Reynolds  nunber  form  of  the  model  other 
inprovements  of  the  standard  model  were  also  incorporated  in  this 
study.  Need  for  inprovsnents  to  the  standard  model  have  been 
recognized  for  saae  time  new  by  severed  researchers  [62,64,69] .  As 
pointed  out  by  Launder  169] ,  the  tern  that  is  most  loosely 
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approximated  in  the  k-e  model  is  e  .  Primarily  because  none  of  the 
terms  in  the  transport  equation  for  e  can  be  measured 
experimentally,  enpirical  approximations  are  difficult.  TO  htprove 
the  modeling  of  the  e  transport  equation  for  accelerated  shear 
flows,  Hanjalic  and  Launder  [62]  "sensitised"  the  equation  to  the 
diffusion  of  energy  effects  by  considering  the  importance  of 
near-wall  turbulence  and  the  influence  of  normal  strains.  The  term 
added  to  the  e differential  equation  is, 
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In  the  above  Ug  is  the  velocity  normal  to  the  streamline,  n  is  the 
co-ordinate  normal  to  the  streamline  and  Rq  is  the  radius  of 
curvature  (Pig.  (HI) ) .  Expressions  for  aug/an  and  P^.,  developed 
in  appendix  (B)  are. 


3ug 

3n 


sin (26) 


f  _  I!i  ] 

l.  3v  3x] 


+ 


cos  (26) 


3u 

3y 


+  •—  ] 
3x  j 


(3. A. 11) 
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(3.4,12) 


In  Bqn.  (3.4.11)  6  is  the  angle  between  the  velocity  vector  that  is 
tangent  to  the  streamline  and  the  x-axis  (Figure  (Bl)}. 

In  the  present  research,  including  the  term  represented 
by  Eqn.  (3.4.10)  resulted  in  some  iitpruveraeur.  in.  the  results  over 
the  standard  model  for  the  recirculating  ?.cne.  However,  it  was 
eventually  discarded  in  favor  of  the  more  comprehensive  approach  of 
including  curvature  effects  by  modifying  ,  as  discussed  in 
Section  (3.4.4).  Results  of  reat tachroent  lengths  and  heat  transfer 
calculated  using  Eqn.  (3.4.10)  in  the  standard  k-e  model  are 
presented  in  a  later  chapter  and  are  compared  with  results  iron 
other  modifications. 

Solution  of  the  equations  of  motion,  emcgy  and 
turbulence  presented  in  this  section  require  specification  of  inlet 
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and  boundary  conditions.  Detailed  discussions  of  these  conditions 
are  presented  in  a  separate  section. 

3.4.4  Incorporation  of  Algebraic  Reynolds  Stress  and  Scalar 

Flux  Closure 

In  this  section  functional  relationships  are  derived  for 
the  coefficients  Cy  and  Pr^  that  are  presented  in  Eqn.  (3.4.6). 
Ibis  procedure  involves  obtaining  transport  equations  for  the 
Reynolds  stresses  and  scalar  fluxes.  These  equations  are  then 
closed  by  the  so  called  algebraic  stress  closure  procedure.  The 
algebraically  closed  equations  are  not  directly  enpicyed  in  this 
study,  but  the  interesting  relationship  between  the  shear  stress 
and  local  rate  of  strain,  as  obtained  from  the  .model,  is  ccrpa:  ed 
with  the  k--e  viscosity  model  of  the  preceding  section  to  obtain  new 
relationships  for  C  and  Pr^.  As  a  result  of  this  application  of 
the  algebraic  Reynolds  stress  and  scalar  flux  closure  scheme,  the 
new  expressions  for  and  Prt  in  Eqn.  (3.4.3),  and  hence  the  k-e 
viscosity  model,  is  now  sensitized  to  account  for  the  effects  of 
streamline  curvatures  and  pressure-strain  (scalar)  interactions 
with  the  inclusion  of  'wall-damping'  contributions.  The  inclusion 
of  these  two  effects  are  crucial  for  obtaining  accurate  solutions. 

The  importance  of  streamline  curvature  effects  in 
turbulent  flow  has  been  demonstrated  in  experimental  studies 
MIG, Ill] .  These  studies  indicate  that  curvature  of  the  mean  flow 
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gives  rise  to  an  appreciable  change  in  the  turbulent  flow 
structure,  -which  then  results  in  a  direct  influence  on  turbulent 
transport  of  heat  and  momentum  [79] .  Attempts  to  include  the 
effects  of  streamline  curvature  in  numerical  calculations  have  in 
the  past  been  based  on  Brashaw's  [68]  modification  of  the  mixing 
length  theory.  However,  because  of  difficulties  in  specifying 
length  scales  and  an  empirical  constant  in  Bradshaw's  equation, 
Leschriner  and  Rodi  [64] ,  Ljuboja  and  Rcdi  [63]  and  Humphrey  and 
Pourhmacii  L£5]  found  inclusion  of  streamline  curvature  effects  in 
the  k-  e  mciel  through  to  be  more  practical.  In  the  present 
study  a  more  generalized  procedure  of  l£f,chziner  and  Rodi  [64j  is 
adopted  for  describing  ,  in  streamline  co-ordinates,  for  the 
recirculating  region.  7he  expression  for  is  generalized  to 
include  'wall-demping1  effects  in  the  pressure-strain  term  and  also 
the  local  equilibrium  of  turbulence  energy  is  not  assumed  (i.e.. 
P/s  »  l.  where  P  is  the  production  of  turbulent  kinetic  energy) . 
In  the  redeveloped  region,  a  generalized  expression  for  is  also 
developed  in  cartesian  coordinates,  similar  to  the  less  rectricted 
procedure  of  Ljuboja  and  Pndi  [63].  In  both  regions,  the  turbulent 
Prandti  number,  Prp,  is  similarly  functionalized.  This  generalized 
procedure  of  improving  the  approximations  of  both.  and  Pr^  in  the 
entire  flow  field  (recirculating  and  redeveloping)  is  unique  to  the 
present  research. 


54 


Ihe  general  outline  of  the  derivation  is  presented  below, 
tedious  algebraic  details  are  presented  in  Appendix  (A). 
Proceeding  from  the  Navier-Stdkes  and  energy  eqiiations ,  the 
equations  that  represent  the  transport  of  Reynolds  stresses,  , 
and  scalar  fluxes  u^0'  ,  respectively,  are: 
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Terms  X,  II/  III/  and  IV  are  often  referred  to  as  the  generation 
(production)/  dissipation,  pressure-strain  (scalar)  and  diffusion 
terms,  respectively. 

For  complete  closure  of  Eqns.  (3.4.13)  and  (3.4.14),  terms  II, 
III  and  IV  need  to  be  approximated.  Launder  [89]  and  Sattiaraweera 
[78]  discuss  in  detail  the  mathematical  procedures  involved  in 
closing  these  unknown  gradient  terms.  This  procedure  is  often 
referred  to  as  the  Reynolds  stresses  (scalar  fluxes)  closure 
scheme.  Bgns.  (3.4.13)  and  (3.4.14)  with  the  closure 
approximations  included  are  then  solved  with  the  transport 
equations  for  mass,  momentum,  energy,  k  and  e  to  obtain  the 
turbulent  flow  field.  A  somewhat  more  restrictive  closure  scheme 
has  been  developed  for  approx iiu&tiiiCi  Eqns.  (3.4.13)  and  (3.4.14) 
and  for  simplifying  the  unknown  gradient  terms  in  Signs.  (3.4.13) 
and  (3.4.14)  151].  In  this  procedure  the  convective  minus  the 

diffusive  transport  terms  are  approximated  by  algebraic  equations 
and  terms  II  and  III  are  approximated  using!  the  same  p-ccedure  as 
that  discussed  by  Launder  [69]  and  Samaraweera  [78] .  It  is  this 
procedure,  often  referred  to  as  the  algebraic  stress  (scalar  flux) 
closure  scheme,  that  is  incorporated  in  the  pie-sent  study  so  as  to 
derive  functional  relationships  for  Cu  and  Prt.  Tht  proo^dure  to 
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(3.4.17) 
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Closure  of  4>j.j  and  Cj_j  results  in  a  set  of  algebraic  equations  for 
the  stresses  ujuj  .  Closure  of  the  pressure-strain  term  is  as 
follows  [693  : 


*ij  "  (1)  +  *ij (2)  +  ^ij (lw)  +  ♦ij (2w) 


(3.4.18) 


where  the  w  subscript  terms  are  the  'wall-damping*  corrections  to 
the  pressure  due  to  turbulence  and  turbulence  and  mean  rate  of 
strain  interactions.  Each  of  the  terms  in  Eqn.  (3.4.18)  are  then 
approximated  as  [17,69] , 
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In  the  above, 

Vi  1JT 

and  C^,  Cs2,  Csiw  and  Cg2w  are  enpirical  constants,  their  values 
are  presented  in  Table  '  (3.1) .  Ihe  near-wall  function,  f (£/»2) » 
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takes  on  a  value  of  unity  near  a  wall  and  zero  remote  from  a  wall. 
Specification  of  f  (  Vxg)  is  as  follows:  For  flews  in  the  presence 
of  a  single  wall  [17] , 


f(£/x2)  «  lQ/y  ,  SLq 
■  3/2 

£(Wx  )  -  ^  >£-_ 
w 
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(3. A. 20) 

where,  is  obtained  by  imposing  the  following  condition, 
f(£/x2)-*  as  y  *  0  {Linear  Sublayer) 

Per  the  'log-law'  Region  (Inertial  Sublayer), 
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(3. A. 21) 


and,  1/2  ' 

•**t  *  **  tp  V  '  (3. A. 22) 

And  fran  the  k-e  viscosity  model. 


p  k2/e  (3. A. 23) 

Subsituting  Dgn.  (3.4.21)  into  Sgn.  (3.4.22),  oonparing  with  Egn. 
(3.4.23)  and  imposing  the  limiting  condition  in  the  viscous 
sublayer,  results  in. 


Using  the  local  equilibruim  assumption,  Cy  »  0.09  [58,1,  the  above 


becomes. 
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1/C  -  0.41 
w 

In  the  present  recirculating  flow  geometries  the  near-wall  function 
should  be  modeled  so  as  to  reflect  the  iiiportance  of  both  side  and 
bottom  walls  in  the  comer  region  of  Fig.  (4}.  Therefore,  in  the 
present  study,  a  superposition  relationship  is  assumed. 
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Rcv-xiting  Bgn.  (3.4.24)  to  correspond  with  the  co-ordinate  systsa 


of  the  present  study.  Fig.  (4),  results  in 
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where,  1  «  0,  for  (D  -  y)  <  d 
x  2  2 

The  only  term  remaining  to  be  closed  is  Eqn.  (3.4.17)  is  the 
dissipation  term.  According  to  Launder,  Reece  and  Rodi  [17] ,  this 
term  is. 


■  I «« 


(3.4.26) 


Subeituting  Bgns.  (3.4.19)  and  (3.4.26)  into  Bgn.  (3.4.17)  results 
in, 
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+  r  [pu -Dij]  +  f  5«  '3-4-27> 

where, 

A’  =  tt-c.j)/[t-i  +  c.l-c.1„f(«/^)] 

and 

A"  =  Cs2w  fW»2)  A,/<1"Cs2) 

Bgn.  (3.4.27)  gives  the  interesting  relationship  between  the  shear 
stress  and  the  local  rate  of  strain;  as  such,  it  is  like  the 
’effective  viscosity'  relationship.  In  the  analysis  to  follow, 
this  concept  is  algebraically  explored  so  # as  to  derive  the 
functional  express  ion  for  Cy  ,  which,  as  pointed  out  earlier,  is 
the  principal  objective  of  this  section. 

A  similar  analysis  to  the  preceding  can  be  performed  for  the 
scalar  flux  transport  equation.  Rewriting  Bqn.  (3.4.14), 

Sj  (^7)  -  F18  +  ♦„  +  D(^TT  (3.4.28) 

where,  Term  I  (production) 

Ae  =  ^18(1)  +  *18(2)  +  ^ie(lw)  +  *i8(2w) 

=  Term  III  (pressure  -  sealer) 

Dtu^Q  )  =  Term  IV  (diffusion) 
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Again,  using  the  closure  procedures  discussed  by  Launder  [69] , 
*10 (1)  "  "C91  k  ui° 

*19(2)'  (3-*'25) 

^iO(lw)  "  'C61w  k  “l9 
♦ieuw)  ■  0 


Values  of  the  empirical  constants,  Cqi,  Ce2»  and  Cqiw  are  given  in 
Table  (3,1), 

Analogous  to  Rodi's  [51]  proposal  for  the  Reynolds  stresses,  the 
convective  and  diffusive  balance  of  scalar  fluxes  can  be  written 


(3.4.30) 


where  Pe  (production  rate  of  6  )  -  -  2  [u^e  31/ Sx^] 


2  "  —  j*1**  —  ■-  ■ 

eg( dissipation  rate  of  e'  )  ■  2  T  (36  /dx^)2 
Substituting  Bgn.  (3.4.30)  into  Bgn.  (3.4.28)  and  making  use  of 
Bgn.  (3.4.29)  results  in. 
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where 
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The  expression  for  Eqn.  (3.4.31)  can  be  sinclified  by  assuming 

local  equilibruim  of  turbulence,  in  which  the  production  and 

.2 

dissipation  rates  of  e  are  assumed  equal,  i.e. , 


V'e  *  1 

Thus, 

*j  *  c«i  *  cei»  ,Wx7)  *  I  <Wc-» 

As  a  consequence  of  the  local  equilibruim  assumption,  the  term 

"2 

containing  c@  and  e’~  in  Bqn.  (3.4.31)  drops  out.  As  such, 

~ 

transport  equations  for  «-6  and  6  are  not  needed.  Bqn.  (3.4.31) 
will  also  be  explored  algebraically  to  yield  a  functional 
relationship  for  the  turbulent  Prandtl  nunber,  Prt,  which  as  stated 
before  is  the  objective  of  this  section. 

In  the  recirculating  region,  adopting  the  procedure  of 
Leschziner  and  Bodi  [64],  Bgns.  (3.4.27)  and  (3.4.31)  are  expressed 
in  strean'ine  coordinates  (s,n) .  As  pointed  out  previously  Figure 
(Bl)  illustrates  the  coordinate  system,  where  coordinate  s  is 
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along  the  streamline  and  n  is  normal  to  the  streamline*  Writing 
these  aquations  for  all  the  combinations  of  stresses  and  scalar 
fluxes  results  in  a  set  of  algebraic  equations.  By  lengthy 
algebraic  elimination  procedures,  unique  expressions  for  the  shear 
stress,  (u'u*  ),  and  the  scalar  flux,  (un6  can  b®  obtained. 

6  U 

Subscripts  's'  and  'n'  refer  to  velocities  tangential  and  normal  to 
streamlines,  respectively*  She  final  expressions  are  presented 
below,  with  the  tedious  algebraic  elimination  details  presented  in 
appendix  (B) . 
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(3.4.33) 


In  the  above  S4  and  are  lengthy  algebraic  express ions  in  terms 
of  turbulent  constants,  iJg/^c,  a*'*  3ug/3n.  Expressions  for  %  and 
gu^Bn  were  given  previously  in  Bqns.  (3.4.12)  and  (3.4.11), 
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i  j.  ->  (a  ^  N, 

1  +  2  K*T  7.  )  R 


(3.4.35) 


The  egressions  for  and  Prt  are  readily  obtained  by  canpairing 
Bgns.  (3.4.32)  and  3.4.33)  with  the  Boussinesq  approximations , 
Bqns.  (3.2.1)  and  (3.2.5)  expressed  in  streamline  co-ordinates. 
The  results  are. 
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where  %2  is  also  a  lengthy  algebraic  expression  given  by. 


Rgn.  (3.4,36)  reduces  to  the  more  restrictive  formulation  of 
Leschziner  and  Rod!  [64]  if  equilibrium  of  turbulence  energy  is 
assumed,  P/e  *  1,  and  if  'wall-dauping1  corrections  in  the 
pressure-strain  terns  are  ignored,  i.e. ,  A  "  *  0. 

In  tiie  preceding  analysis  (and  Prt)  are  modeled  in 

streamline  coordinates  to  reflect  the  sharp  bending  of  the 
streamlines  that  dominate  in  the  recirculating  region.  In  the 
present  study  the  expressions  for  and  Prt,  Bgns.  (3.4.3G)  and 
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(3.4.37),  are  used  with  the  standard  k-e  model  for  predicting  the 
flow  field  in  the  recirculating  region.  In  the  redeveloping 
region,  where  tlie  bending  of  the  streamlines  are  not  so  severe, 
expressions  for  (to  account  for  curvature  effects)  and  Prt  are 
obtained  in  (x,y)  coordinate  system.  The  matheTiatical  procedures 
involved  in  obtaining  these  expressions  for  the  redeveloping  region 
are  analogous  to  the  previous  derivation.  They  however,  result  in 
even  more  tedious  algebraic  detail  and  lengthy  algebraic 
expressions-  Consequently,  the  mathematical  details  are  presented 
in  Appendix  (B).  The  final  expressions  are: 


in  which. 


F  +  BC  +  ~  lcB(l  +  BC) 
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A,B,C/C  ,E,F,  a  and  6  are  lengthy  algebraic  expressions  involving 
turbulent  constants  and  mean  gradients  arvd  are  presented  in 
Appendix  (B) .  The  above  expression  for  ,  Bgn.  (3.4.38),  reduces 
to  the  expression  of  Ljuboja  and  Podi  [63]  if  the 
convective/diffusive  transport  of  u'v*  and  longitudinal  gradients, 
3v/9y/3<i/3y  ,  are  ignored. 

The  expressions  given  above  for  C  ^  ard  Prt,  valid  for  the 
redeveloping  region,  are  used  with  the  low-Reynolds  number  version 
of  the  k-e.  model.  When  using  the  expression  for  ,the  near-wall 
function,  f^  (Dgn.  3.4.6),  in  tne  lcw-Revnolds  number  model  is 
emitted.  This  is  essentially  because  <1  is  already  modeled  to 
reflect  the  near-wall  effects.  The  physical  procedures  involved  in 
dividing  the  flow  field  beyond  rearward  steps  and  pipe  expansions 
into  two  distinct  regions,  for  application  of  either  the  standard 
or  lc*r= Reynolds  nnsbcr  model,  are  described  in  the  next  section. 
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3.5  Two-Pass  Procedure 

Initial  calculations  made  by  simply  employing  the  standard  k-  c 
model  formulation  and  the  corresponding  wall-function  boundary 
conditions  yielded  unsatisfactory  results  for  the  heat  transfer 
coeff icients  in  the  redeveloping  region  [106] .  The  predicted  heat 
transfer  coefficients  in  this  region  were  significantly  higher  than 
the  data.  Using  the  low-Reynolds--nurufcer  formulation  results  in  a 
large  and  erroneous  reduction  in  heat  transfer  coefficients  in  the 
recirculating  region,  but  better  agreement  is  obtained  with  the 
data  in  the  redeveloped  region.  This  behavior  is  understandable, 
because  the  extra  teems  appearing  in  the  low-Reynolds-number  model 
were  a  result  of  tests  made  on  she.,  r  layers  [60] .  Thus,  using  the 

4» 

low- Reynolds  nunber  model  in  recirculating  flews,  upstream  of 
reattachment  will  not  be:  totally  consistent. 

Improvement  of  solutions  for  heat  transfer  in  the  entire  flow 


X.6V4UXJ.CU  a  of 


Tp  -p  1  wf* 


calculation,  the  reattachtent  point  .is  estimated  using  the  improved 
standard  k-  e  model.  In  general,  as  will  be  shown  in  a  later 
chapter,  this  initial  solution  tends  to  deviate  somewhat  from  the 
experimental  results,  particularly  in  the  redeveloped  region,  but 
it  does  predict  fairly  accurately,  with  the  appropriate 
modification  to  C  ,  the  reattachraent  point.  Once  the  reattachment 
point  is  determined,  the  calculation  is  repeated  with  the  improved 
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standard  k-  e  model  employed  upstream  of  reattachment,  and  the 
improved  lcw-Reynolds  number  version  of  the  model  applied 
downstream  of  reattachment.  Figure  (5)  shows  the  solution  domains 
for  each  model  during  the  second  calculation.  In  implementing  this 
procedure  during  the  second  pass,  the  low- Reynolds  number  terms  in 
Bqns.  (3.4.6)  -  (3.4.8)  are  simply  cmittted  upstream  of 
reattachment.  The  boundary  conditions  for  this  procedure  are 
discussed  in  the  next  chapter. 
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3.6  Treatment  of  Initial  Conditions  for  Ccmuts-d  Flews 

The  upstream  conditions  for  each  test  case  coBputed  i:  this 
study  are  summarized  in  Table  3.2.  Unfortunately,  the  wtk  cZ  Awg 
and  Goldstein  IS]  does  not.  include  complete  ute^suraa*nts  of 
upstream  (approach)  profiles  since  their  focus  was  on  hear- transfer 
coefficients.  They  do,  however,  include  measurements  of  approach 
temperature  boundary-layer  profiles.  Their  experiment  was 
performed  utilizing  a  tripped  approach  boundary-layer  with  a  fairly 
short,  but  not  entirely  insignificant,  thermal  entry  length.  The 
approxlTvatioiiS  of  the  inlet  profiles  for  u,  k  and  e  are  described 
below.  A  similar  procedure  as  that  used  to  obtain  the  velocity 
profiles  was  used  to  determine  "the  temperature  profiles,  which  when 
conpared  with  the  data  revealed  no  inconsistencies.  Thus,  serving 
as  a  check  an  the  validity  of  the  procedure  for  the  velocity 
profiles.  The  simulation  procedure  is  presented  below,  with 
algebraic  details  presented  in  Appendix  (C).  The  various  lerxiuj 
that  appear  In  the  equations  below  (e.g. ,  X]  and  i,  )  are 
schematically  presented  in  Figure  (Cl).  Frau  turbulent  flat  piax* 
boundtiry-.layer  theory  1.103]  the  momentum  boundary" layer  thickness, 
at  the  step  is. 


The  corresponding  unheated  starting  length  thermal  boundary  layer 
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thickness  at  the  step  is  [103] , 

W  “  W 

(3.6.2) 

in  which  5  ■  xj"  _  xy",  x^"  and  X2“  an  defined  below.  Qgns. 
(3.6.1)  anc  (3.6.2)  require  that  the.  effective  starting  length  x^ 
be  known.  The  procedure  to  determine  x  _  is  as  follows.  (See  Figure 
(Cl)  for  the  schematics.)  The  momentum  boundary  layer  up  to  the 
trip  is  assumed  laminar,  given  the  starting  length,  as  specified  by 
the  expe  riroa;..  [8],  thus  boundary  layer  thickness  at  the  step  is 
readily  determined.  Assuming  that  imru-diately  at  the  step  the 
bounda-.y  l«s/*sr  is  turbulent,  anc.  equating  the  obtained  laminar 
value  to  the  turbulent  growth  law.  results  in  &r»  effective  starting 
length  for  a  turbulent  boundary  layer  at  the  trip,  xj".  Thus  the 
total  effective  starting  1  ength  ,xy,  is  sum  of  xj/  and  the  distance 
between  the  tr  p  and  the  3tep,  xs,  A  similar  procedure  is  applied 
to  obtain  the  equivalent  tutoulent  thermal  starting  l&gth  at  the 
trip,  %2\  Wuan  th.&  distance  is  used  in  Eqn.(3.6.2),  the  thermal 
boundary  layer  thick iM&s  ccupcres  favorably  with  the  experiments* 
values  [8] . 

Once  the  momentum  and  th.raal  boundary  layer  thick,  are 


deberminod,  the  velocity  and  tssiperatur:?  profiler,  at  the  sfeyp  can 
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# 


be  determined  by  using  the  Abbott  and  Kline  L10J  fit  for  a  tripped 
boundary  layer. 


w 


1/5.6 


(3.6.3) 


Rewriting  the  above  for  the  coordinate  system  used  in  this  study. 
Figure  (4), 
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(3.6.4) 


Ass  lining  that  tenperature  develops  similarly  as  the  velocity, 
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(3.6.5) 


In  the  abo/e  U0  and  T0  are  the  free  stream  velocity  and 
teiperature,  respectively,  and  Tw  is  the  wall  taiperature.  Inlet 
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profiles  for  k  and  e  are  obtained  by  a  piecewise  linear  fit  of  the 
flat  plate  boundary  layer  data  of  Klebanoff  1112]  as  presented  in 
Hinze  [16].  In  the  free  stream,  k  is  assigned  a  value  of  0.1 
percent  of  the  free  stream  velocity  (i.e.,  k®0.001U^).  In  initial 
calculations  [106]  it  was  discovered  that  the  level  of  turbulence 
intensity  in  the  mixing  region  of  the  fit  led  to  unreasonable 
reattachment  lengths  when  coopered  with  the  data.  lb  eliminate 
this  difficulty  it  was  necessary  to  increase  the  level  of  turbulent 
kinetic  energy  in  this  region  approximately  three-fold.  This 
higher  level  of  turbulent  energy  i/>  not  completely  arbitrary,  since 
the  presence  of  the  trip  in  the  experiments  would  seem  likely  to 
produce  turbulence  levels  persisting  downstream  of  the  trip  higher 
than  those  of  the  corresponding  Klebanoff  data  for  the  same 
Reynolds  number.  Subsequent  calculations  studied  the  influences  of 
the  trip  (which  is  situated  approximately  11  trip-heights  upstream 
of  the  step)  on  approach  conditions  at  the  step  and  on  the  net 
effect  on  reattachroent  lengths.  The  results  indicate  that  in  order 
to  have  reattachment  length  similar  to  the  experiment,  the 
turbulent  intensity  at  the  step  xust  be  equivalent  to  the  intensity 
found  at  11  step-heights  downstrean  of  a  step,  which  is 
significantly  much  more  enhanced  than  the  Klebanoff  [112] 
intensities.  Results  of  these  studies  are  presented  in  a  later 
chapter. 
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In  the  data  of  Klebaroff  [112j  the  turbulent  kinetic  energy 
(and  c  )  are  normalized  by  relations  containing  the  friction 
velocity,  .  For  a  turbulent  flat  plate  uT  is, 


nl  a  0.0066  U2  Ra“1/6. 
r  o  GUO 


(3.6.6) 


The  heat  transfer  results  of  Aung  and  Goldstein  [8]  are 
normalized  by  the  Nusselt  number  just  upstream  of  the  step  (Nutj) . 
This  quantity  is  not  provided,  and  so  must  be  calculated.  This 
quantity  is  readily  obtained  from  turblent  flat  plate  theory  1103] , 


0.8  ! 

S/10l  -1/9 

0.0287  Reg  ?r  J 

_ l  -  (5/*)  J 
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(3.6.7) 


0.169  [13.2  Pr-10.16]  +  0.9 

Because  the  trip  is  at  a  distance  of  IX  trip-heights  upstream  of 

the  step,  it  is  necessary  to  enhance  the  result  of  Bgn.  (3.6.7)  by 

a  factor  of  about  1.4  [8] .  The  reason,  as  will  be  shown 

graphically  in  a  later  chapter,  is  that  the  heat  transfer 


coefficient  beyond  a  step  (in  this  case  the  trip  is  the  step)  never 
approaches  the  exact  flat  plate  value. 

Fbr  the  back-step  experiments  of  ttoss  and  Baker  [9]  and  Eaton 
and  Johnston  [20J  all  approach  profiles  are  reported,  except  for 
the  turbulence  dissipation,  which  is  obtained  from  the  relationship 


below. 
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3/2 

e  „  -  £  -  0.005  d/2 

A  O 


(3.6.8) 


Results  are  not  very  sensitive  to  the  choice  of  l  q. 

Ffor  the  pipe  expansions  of  Zemanick  and  Dougall  [39 J  and 
Sparrow  and  O'Brien  [41] ,  fully  developed  turbulent  inlet  profiles 
are  specified  for  the  velocity  and  temperature.  These  are. 
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(3.6.9) 


f'X  c  i  n  \ 


•The  profiles  for  k  and  e  are  obtained  by  a  piecewise  linear  fit  of 
the  data  of  Laufer  [113].  Corresponding  to  Bgn.  (3.6.6)  for  the 
flat  plate,  the  friction  velocity  for  fully  developed  turbulent 
pipe  flow  is, 

2  _1M 

u  -  0,04  D  Sa_  (3.6.11) 


The  heat  transfer  coefficients  for  the  Zemanick  and  Dougall 
experiment  [39]  require  calculations  of  the  bulk  temperature,  it 
can  be  easily  verified  that  the  inlet  bulk  temperature  is. 


t  „  ZuTdA 
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T  -  0.833  (T-T) 
w  w  » 


(3.6.12) 
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Presented  in  this  study  are  correlations  of  back-step  heat 
transfer  coefficients  for  various  Reynolds  number.  TO  rule  out  the 
effect  of  approach  boundary  layer  thicknesses  on  heat  transfer, 
constant  neroentum  boundary  layer  thickness  is  assuioed  for  ail 
Reynolds  numbers.  (In  the  present  study  the  effect  of  approach 
boundary  layer  thickness  on  reattachment  length  is  shewn  to  be 
negligible.)  Also,  the  thermal  and  momentum  boundary  layers  are 
assumed  tc  have  the  same  starting  length.  The  momentum  boundary 
layer  thickness  at  the  step  is  assumed  to  be, 

Su  »  C  H,  where  C  is  a  constant  chosen  as  described  beiew 

(3.6.13) 


Fran  turbulent  flat  plate  theory  [103] ,  the  corresponding  starting 


length  is, 

° 


•J..25 

2.6219  <Su(pUQ/u) 


(3.6.14) 


Using  Eqri.  (3.6.2),  with  £  »  0,  5i(xl)  at  the  step  cart  be  obtained. 
In  Bgn.  (3.6.13)  the  values  of  C  examined  are  0.4  <  C<  2.0.  For 
the  correlations  C  is-  taken  to  be  1.5.  Discussion  on  the  effects 
of  the  range  of  C  on  heat  transfer  coefficients  is  presented  in  a 
later  chapter. 
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3.7  Final  Remarks 

In  this  chapter  closure  schemes  for  the  turbulent  viscosity 
and  thermal  diffusivity  terms  in  Bgns.  (3.3.2)  -  (3.3.4)  have  been 
presented  together  with  some  theoretical  items  relevant  to  the 
computational  process.  In  the  present  study  of  turbulent 
recirculating  heat  transfer,  approximations  of  the  viscosity  by  a 
single  mixing  length  scale  proved  to  be  inadequate.  Variations  and 
improvements  of  the  k-e  model  were  then  tested.  Hie  final 
calculations,  the  results  of  which  are  to  be  presented  in  a  later 
chapter,  employ  the  two-pass  procedure  with  the  improved  version  of 
the  standard  k-e  model  dismissed  and  the  low-Reynolds  number  form 

V>J,  UiC  HKAU • 

Direct  application  of  the  relatively  complicated 
"second-order" ,  Reynolds  stresses  and  scalar  fluxes  closure  schemes 
have  not  been  attempted  in  this  study.  It  was  not  necessary 
because  the  two-pass  method,  as  seen  In  a  later  chapter,  adequately 
predicts  heat  transfer  coefficients  for  engineering  purposes. 
Also,  using  the  Reynolds  stresses  and  scalar  fluxes  procedure  would 
require  initial  conditions  for  each  of  the  stresses  and  scalar 
fluxes.  There  is  a  lack  of  measurements  of  these  correlations  for 
approach  boundary  layers.  Thus,  the  accuracy  of  the  model  cannot 
be  thoroughly  tested  in  the  absence  of  well  defined  initial 
conditions.  Moreover,  engineering  applications  favor  simple 
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modeling  procedures  and  the  present  two-pass  k-e  calculation  is  a 
simpler  approach  than  carplete  Reynolds  stress  closure  modeling. 

In  the  next  chapter  the  numerical  procedures  used  to  solve  the 
set  of  differential  equations,  Eqns.  (3.3.1)  -  (3.3.4)  and  Eqns. 
(3.4.7)  ••  (3.4.8),  are  presented. 


CHAPTER  IV 


COMPUTATIONAL  PROCEDURES 

4.1  General  Remarks 

The  computational  procedure  outlined  herein  is  an  outgrowth  of 
the  TEACH-T  program  described  by  Gosraan  [101] .  The  program  in  its 
original  form  features  a  variable,  staggered  grid  arrangement  and 
.solves  the  hydrodynamic  equations  described  in  Chapter  III  by 
finite  differencing  procedures  based  on  flux  balances.  In  the 
present  study,  a  temperature-field  solution  routine  together  with 
improvements  upon  differencing  of  the  convective  terms  were 
incorporated  into  the  code.  Hie  flow  considered  is  two-dimensional 
and  the  equations  are  written  in  generalized  cylindrical 
co-ordinates  (x,r)  and  can  be  readily  converted  into  Cartesian 
co-ordinates  (x,y)  by  changing  an  index  in  the  program. 


8'J 

4.2  Numerical  Method 

In  this  present  work,  attention  is  focused  on  solving  the  full 
ITavier-Stokes  and  energy  equations  for  steady  two-dimensional 
flews.  The  numerical  procedure  used  was  developed  by  Gosnan  and 
coworkers  LlOj  and  is  an  outgrowth  of  the  SIMPLE  procedure  of 
Patankar  [95J  developed  for  boundary  layer  flows.  Details  of  the 
finite  differencing  procedures  are  discussed  very  thoroughly  by 
Patankar  [95,96];  therefore,  only  a  brief  overview  of  the  solution 
procedure  is  presented  herein. 

4.2.1  TEACH-T  Solution  Algorithm  with  Heat  Transfer 

Ihe  equations  in  Chapter  III,  Bqns.  (3.3.2)  -  (3.3.4)  and 
Bgns.  (3.4.7)  -  (3.4.8),  can  be  written  generally  as, 

4(prH*  ■ 


=  0 


where 


g  =  u,  v,  T,  k  or  e 


(4.2.1) 


and 


r  = 


yeff  or  reff 


In  Cartesian  coordinates,  r  =  1  and  ^/5r  =  3/3y.  is  the  source 
term  and  represents  the  mechanism  for  the  generation  or  destruction 
of  <f>.  However,  also  represents  all  other  terns  that  cannot  be 
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conveniently  expressed  by  the  convection  and  diffusion  terms  of 
Eqn.  (4.2.1),  for  example,  the  pressure  gradient  term  (3p/3x  )  in 
the  momentum  equations.  The  mass  flow  rate  in  Eqn.  (4.2.1)  must 
satisfy  the  continuity  equation.  This  condition  is  satisfied  by 
updating  the  pressure  solution  so  that  the  resulting  velocity  field 
satisfies  the  continuity  equation.  The  solution  for  the  pressure 
field  is  discussed  later. 

The  discretized  equations  are  obtained  by  integrating 
Dgn.  (4.2.1)  over  a  control  volume.  Figure  (6)  illustrates  the 
control  volumes  and  grid  points  used  in  the  TEACH-T  program.  'The 
control  volumes  are  indicated  by  the  broken  lines,  and  the  shaded 
area  represents  cne  such  control  volume.  The  grid  points  are 
placed  at  the  geometric  centers  of  the  control  volumes,  the  solid 
lines  represent  the  grid  lines.  This  procedure,  often  referred  to 
as  the  staggered  grid  arrangement,  aids  in  numerical  stability 
1101]*  Figure  (6),  for  simplicity,  shows  grid  lines  of  constant 
spacing:  however,  in  this  study  non-uniform  grid  distributions  are 
used.  One  such  distribution  is  shewn  in  Figure  (7).  (It  should  be 
pointed  out  that  in  Figure  (7)  the  magnifications  .of  the  vertical 
and  horizontal  scales  are  not  of  the  soma  order. )  This  grid 
distribution  is  devised  to  concentrate  grid  lines  near  the  planes 
of  the  walls.  As  will  be  discussed  later  various  near-wall  grid 
densities  were  tested  to  obtain  grid- independent  heat  transfer 
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results.  Figure  (6)  also  indicates  the  treatment  of  the  near-wall 
control  volumes;  for  instance,  at  a  wall  grid  points  are  placed  at 
the  centers  of  the  control  volume  faces.  The  advantage  of  this 
procedure  is  that  it  allows  the  treatment  of  different  boundary- 
conditions.  For  instance,  the  value  of  4>  at  the  boundary  or  the 
flux  across  the  control  surface  can  be  conveniently  specified. 

With  the  control  volumes  defined,  Eqn.  (4.2.2)  can  be 
integrated  for  each  control  surface.  This  procedure  is  exemplified 
by  considering  an  isolated  node  point,  P,  as  shewn  in  Figure  (8). 
Eqn.  (4.2.1)  then  becomes, 

/  [(eru*  -  rr  ||  ^  -  (  Pru<f>  -  rF  ||  )  J  dr 

+  J'l  (  Prv<*>  -  rr  ||-  )  -  (prv<J>  -  rr  ||  )  ]  dx 

wL  n  s  J 

-  f  S .  r  dr  dx  «  0  (4.2,2) 

V  9 


By  having  'J'  represent  the  total  convection  and  diffusion  flux  in 
a  given  direction,  and  evaluating  the  integrals  by  the  mean  value 
theorem,  Bgn.  (4.2.2)  reduces  to, 

J  A.  -  J  A  +  J  A,  -  J  A  +  S,  A  V  «  0 

ww  ee  ss  nn  <? 


(4.2.3) 
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Approximation  of  the  fluxes,  Ji,  in  the  original  TEACH-T  code  is 
accomplished  by  the  hybrid  differencing  schem  [95] .  In  this  study, 
the  hybrid  differencing  scheme  is  replaced  by  the  quadratic 
differencing  schane  [98] .  Both  of  these  differencing  schemes, 
together  with  the  justification  for  the  change,  will  be  discussed 
in  the  next  two  sections.  Often  the  source  term  contains  the 
variable  <J>,  and  can  be  linearized, 

S ^  =  Sc  +  Sp  <J>p  (4.2.4) 

where,  Sc  is  the  part  of  that  has  no  <J>  dependency  and  Sp  is  the 
coefficient  of  $p.  with  these  approximations  Bgn.  (4.2.3)  can  be 
discretized  [95] . 

The  final  discretized  equations  are  only  qussi-1  inear ; 
for  instance,  the  coefficients  may  depend  on  values  of  4>  and  also 
on  other  dependent  variables.  As  a  result,  the  final  solution  is 
obtained  iteratively.  Convergence  can  be  attained  by  a  series  of 
iterations.  However,  in  most  cases,  convergence  is  not  simply 
obtained  as  a  result  of  successive  iterations;  for  instance,  the 
values  of  4>  may  oscillate  or  diverge.  Controlling  the  changes  in  <J> 
is  accomplished  by  the  under-relaxation  technique  [95] .  Fbr 
example, 

4*P  =*<£?+  a  [changes  in  <f>p  in  the  current  iteration] 

* 

where  p  is  the  value  of  p  fran  the  previous  iteration.  To  reduce 
the  change  in  $p,a  ranges  fran  0  to  1.  Convergence  is  achieved 


B4 


when 

What  finally  remains,  is  a  solution  procedure  for  the 
resulting  set  of  quasi-linear  algebraic  equations.  h 

one-dimensional  linear  discretized  set  of  equations  can  be  solved 
efficiently  by  the  tridiagonal  matrix  algorithm  f  95] .  For 

two-d imens ional  problems,  a  line-by-line  solution  technique  195 J 
incorporating  the  tridiagonal  algorithm  is  utilized  in  the  TEACR-T 
code.  In  this  procedure,  the  equations  for  all  the  <p  values  along 
one  grid  line  are  considered,  and  the  neighboring  line  values  of  <J> 
are  subsituted  from  the  best  known  estimate.  One  iteration  of  the 
line-by-line  technique  is  completed  when  all  the  lines  in  a  chosen 
direction  are  swept. 

The  pressure  is  updated  by  using  the  equation  of 

continuity  in  conjunction  with  a  ’Poisson-like'  equation  195] .  For 
instance,  the  velocity  is  expressed  in  terras  of  pressure  (for 

simplicity  looking  at  only  the  w-direction  of  Figure  (3)), 


9u , 


w 


“'.*3757^7  ■  V 


r  a  '»  c  > 


in  which  p  is  pressure  correction  and  the  term  /<KPW-Pp)  is 
obtained  from  the  solution  of  the  u-manentum  equation.  The 
integrated  form  of  the  continuity  equation  is, 

(puA)e  -  (puA)w  +  (pvA)n  -  (pvA)s  *  0 


(4.2.6) 
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Substituting  for  the  velocities  in  Eqn.  (4.2.6)  by  using 
relationships  like  that  of  Eqn.  (4.2.5)  for  ue,  uw,  un  and  us 
result  in  a  discretized  equation  for  the  pressure  correction.  The 
pressure  obtained  frcm  solution  of  this  equation  is  then  used  to 
correct  the  velocities  to  obtain  velocity  corrections,  u  and  v. 

The  numerical  solution  procedure  for  the  TEAOi-T  code  is 
as  follcvs: 

1.  Initialize  fields  for  all  variables. 

2.  Assemble  coefficients  of  momentum  and  energy 

equations  and  solve  for  new  variables  (u*fv*,  etc.) 
by  using  previous  pressure. 

3.  Solve  for  p  in  pressure  correction  equation. 

4.  Update  all  fluid  flow  variables,  i.e., 
p  =  p*  +  p,  u  =  u*  4-  u,  etc. 

5.  Assemble  coefficients  and  solve  for  all  other 

values,  i.e. ,  k,  and  c. 

6.  Test  for  convergence. 

7.  If  convergence  not  attained,  continue  solution 

procedure  frcm  step  (2)  until  convergence  is 

achieved. 

In  this  study  convergence  of  the  solution  is  considered 
satisfactory  when  the  normalized  residuals,  R,  of  each  <f>  equation, 
summed  over  the  entire  calculation  domain,  are  smaller  than  0.005. 
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Dor  example,  for  4>=u, 


lu  -  momentum  equation 

R  *  —j -  (4.2.7) 

j  p  uinlet  d  A 

In  other  words,  convergence  is  achieved  if  Bgn.  (4.2.2)  is 
satisfied.  In  the  TEACH-T  code  there  are  three  basic  termination 
procedures;  first,  if  after  20  iterations  and  the  maximum  residual 
source  is  excessively  large,  then  it  will  be  impossible  bo  have 
converged  solutions,  second,  if  convergence  is  obtained  and  third, 
if  the  nunber  of  iterations  has  reached  the  maximum  value 
specified.  In  this  study,  the  maximum  number  of  iterations  has 
been  set  at  400. 

As  mentioned  previously,  the  TEACH-T  computer  code  in  its 
original  form  did  not  include  temperature  calculations.  However, 
solution  of  the  energy  equation,  Bqn.  (3.3.4)  is  easily 
incorporated  in  the  code  by  mimicking  the  numerical  procedures  of 
say  the  u  -  momentum  equation.  To  validate  the  code  initial 
calculations  ware  performed  for  laminar  heat  transfer  over  a  flat 
plate  and  more  recently,  in  as  yet  unreported  work  by  the 
Department  of  Mechanical.  Engineering,  Howard  University, 
laminar  heat  transfer  flows  over  a  cavity. 


for 
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4.2.2  'IEACH-T  Hybrid  Differencing  Scheme 

The  hybrid  differencing  procedure  used  in  the  solution 
algorithm  is  discussed  in  detail  by  Patankar  [95] .  Ihis  procedure 
uses  both  the  central  and  upwind  differencing  schemes. 
Approximations  for  the  convective  terms  in  the  transport  equations 
are  as  follcws:  for  face-Peclet  number,  Pe,  (pu^x^/f,  convective 
to  diffusive  transport)  less  than  2,  central  differencing 
approximations  are  used  and  for  face-Peclet  nunber  greater  than  2, 
upwind  differencing  approximations  are  used.  Hie  diffusion  terms 
are  calculated  by  central  differencing  approximations. 

From  Figure  (8),  central  differencing  of  $  e  is  the 
average  of  and  <j>p,  and  upwind  differencing  uses, 

4>e  =  4>p  ,  for  u  >  0 
4>  *  <pK  ,  for  u  <  0 

By  using  the  above  approximations  differenced  equations  can  be 
obtained,  as  presented  below.  Centred  differencing  of  the 
diffusion  term  is, 


Sx2 


2  d>.  +  <£.  , 

Yi  ^i-l 


(Ax)J 


(4.2.8) 
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And  the  convection  term  is. 
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u  <  0,  Pe  >  2 


(4.2.9) 


(4.2.10) 


(4.2.11) 


Leonard  L98J  showed  that  for  convection  dominated  flows  the 

hvrHr*i  H  nnrvfiH^Trp  ic  1  \  ?  mavrrinal  rnr^vi  /->*-»  ■*-*» 

—  —  —  *-  —  —  WMr“k*v*  wv  U^TTXIIU 

differencing.  He  further  demonstrated  that  the  hybrid  procedure 
severly  restricts  the  grid  densities  (spacings),  consequently, 
limiting  the  numerical  accuracy,  and  also  that  the  upwind 
differencing  scheme  is  a  main  source  of  numerical  diffusion.  To 
eliminate  same  of  these  discrepancies  Leonard  [90J  proposed  the 
quadratic  differencing  scheme. 

4.2.3  New  Quadratic  Differencing  Scheme 

The  quadratic  scheme  of  Leonard  198J  seeks  to  ccsnbine  the 
accuracy  of  the  central  difference  scheme  with  the  stability 
property  of  the  upwind  scheme.  Rather  than  using  the  linear 
interpolation  scheme  of  the  central  difference  procedure,  a 


parabolic  polynomial  interpolation  scheme  is  employed.  The  control 
volume  surface  values  are  obtained  by  fitting  a  parabola  to  the 
values  of  <P  at  three  consecutive  nodal  positions,  two  located  on 
either  side  of  the  face  considered  and  the  third  located  at  the 
next  node  in  the  upstream  direction.  As  such,  consideration  of  all 
the  faces  for  a  point  P,  as  in  Figure  (8),  would  require  a  9-point 
star  scheme  in  which  the  sign  of  all  the  face  velocities  must  be 
taken  into  consideration.  The  quadratic  approximations  for  say 
face  'e'  (Figure  (8) )  are  as  follows  [98 J : 


♦.  *  7  »i  *  +i .  i>  -  r  «i- 1  *  *  i  - 2  V-  V0  (4-2-12) 

♦e  -7  »i  *  ♦,*!>  -  ?  »i  *  -  2  ♦i.p.  u.<0  (4-2-13> 


Similar  expressions  can  be  written  for  the  other  three  faces.  In 
Eqns.  (4.2.12)  and  (4.2.13),  the  first  bracketed  terms  represent 
the  linear  centred  differencing  approximations,  and  the  second 
bracketed  terms  represent  the  correction  of  <p  by  upstream  weighted 


.a.aa;  ana  inco 


the  TEACH-T  code  (described  previously)  led  to  an  increase  in  the 
number  of  iterations  needed  for  convergence  try  as  much  as  100 
percent.  This  is  probably  due  to  the  loss  of  diagonal  dominance  as 
a  result  of  the  extra  terms  in  the  quadratic  differencing 
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was  only  applied  to  the  momentum  and  energy  transport  equations. 
This  was  because  the  recent  study  of  Leschziner  and  Rodi  [63]  has 
shown  that  solutions  of  k  and  e  are  insensitive  to  the 
discretization  scheme,  particularly  in  the  high  shear  regions 
bordering  recirculating  zones.  However,  in  order  to  obtain  grid 
independent  solutions,  grid  refinaner.ts  were  necessary  in  the 
near-wall  region;  for  instance,  about  a  10  percent  increase  over 
what  was  needed  when  all  the  transport  equations  incorporate  the 
hybrid  difference  scheme.  Recently,  by  extensive  numerical 
experimentations,  Han,  Humphrey  and  Launder  [100]  showed  that  for 
grid  distributions  less  than  (22  x  22)  in  cavity  driven  flows, 
quadratic  differencing  approximations,  for  the  convective'  terms 
should  be  used  for  all  variables.  As  a  result,  subsequent 
calculations  in  this  stix3y  did  incorporate  quadratic  differencing 
in  all  of  the  transport  equations.  With  this  procedure  grid 
independent  solutions  were  very  easily  obtained  with  about  a  (25  x 
28)  grid  distribution. 

As  a  final  remark,  seme  of  the  advantages  gained  by  using 
the  quadratic  differencing  procedure  have  been  demonstrated  by 
Leonard  [98] .  !he  list  be  lew  represents  a  brief  summary  of 
Leonard's  [98]  study,  where  quadratic  differencing  of  the 
convective  terms  results  in: 
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U)  Reduction  of  the  presence  of  false  diffusion, 

(ii)  More  stability  at  high  Peclet  numbers. 

(ni)  Less  grid-dependent  solutions  than  the  upwind 
scheme. 

As  a  result  of  the  above  conclusions  this  improved  scheme  is  much 
more  realiable  in  testing  the  validity  of  turbulent  models. 


4.3  Treatment  of  Boundary  Conditions 


In  order  to  avoid  an  excessively  fine  grid  at  a  wall, 
functional  relationships  are  specified  to  describe  the  rapidly 
varying  quantities  in  the  wall  region.  The  procedure  of  obtaining 
these  functions  are  two  fold:  the  values  for  sane  variables  (e.g. , 
u,v,T)  at  the  nodes  adjacent  to  the  wall  (near-wall  nodes)  are 
directly  obtained  from  'wall-function'  relations,  and  where 
'wall-function'  approximations  are  not  strictly  valid,  approximate 
solutioas  of  the  transport  equations  for  the  other  variables  (e.g. , 
k  and  e  )  are  solved  in  the  near-wall  region.  In  this  study,  as 
mentioned  previously,  two  forms  of  the  k-e  model  are  employed  in  a 
two-pass  procedure.  The  standard  model  is  applied  to  the  entire 
flow  field  in  a  first  pass  procedure  to  establish  the  recirculating 
zone,  and  in  a  second  pass  the  standard  model  is  maintained  up  to 
reattachment  and  the  lew- Reynolds  number  version  is  applied 
downstream  of  reattachment.  The  boundary  conditions  for  each  model 
require  slightly  different  treatments,  consequently  in  the 
following  sub-sections  each  procedure  is  analyzed  separately. 

4.3.1  Single  Pass  Procedure 

The  standard  k- e  model,  is  discussed  in  Chapter  III,  is 
valid  for  regions  outside  the  buffer  zone  (Figure  (C2))  where  the 
flow  is  fully  turbulent.  Beyond  this  region,  13  <  y* (uTy/v)  <250, 
is  the  'log-law  region'  (Figure  (C2)),  and  in  this  region 


functional  relations  are  assigned  to  the  velocities  and  temperature 
158].  These  relationships,  often  referred  to  as  ’wall-functions', 
are  approximate  solutions  of  the  transport  equations.  The 
functioned  relationships  for  the  entire  near-wall  region  (Figure 
(C2) )  for  momentum  and  energy  transfer  are  obtained  by  integrating 
the  governing  boundary  layer  transport  equations  in  each  region 
with  the  appropriate  assumptions  applied.  The  procedure  is  briefly 
described  below  and  in  Appendix  (C)  the  mathematical  procedure  is 
carried  out  for  the  u-mcmentum  equation.  In  the  linear  sublayer, 
<  5,  the  flow  is  assumed  laminar  and  the  solution  of  the 
Navier-Stckes  and  energy  equation  is  obtained  by  assuming  the 
shear-stress  (and  heat  flux)  is  constant  throughout  this  region. 


(4.3.1) 

(4.3.2) 


where 

Tt  E  VH  UT 

In  the  buffer  region,  5  <yf  <  13,  both  laminar  and  turbulent  shear 
stresses  (and  scalar  fluxes)  are  important,  so  approximate  solution 


y5 

of  the  wall  layer  transport  equations  results  in  an  integral 
equation.  In  the  log-wall  region,  13  <  y*  <  230,  the  flew  is  assumed 
fully  turbulent,  and  the  Boussinesq  approximations  of  Chapter  III 
are  applied.  Solution  of  the  transport  equations  result  in, 

ui+  *  k  Zn  xi  +  C*  K  *  °*42  (4.3.3) 

T+  =■  ~  in  x*  +  D  (4.3.4) 

K  1  '■ 


In  the  TEACH-T  code  the  boundary  conditions  are  expressed  through 
wall  shear  stress  and  heat  flux.  Bor  the  standard  k-e  jiodel  the 
'wall-function'  relations  are  summarized  below,  and  they  are  those 
proposed  toy  Launder  and  Spalding  158] .  These  relationships  are 
modifications  of  the  boundary  layer  results  (presented  above)  so  as 
to  reflect  the  presence  of  recirculating  flew. 

{ i)  Mean  velocity:  the  vail  shear  stress, 

_1 

T»  ■  »„  “p  c-u  kp  —T— r  •  E  ■  s-8  <4-5-« 

M  Jin  (Eyp  ) 

where  subscripts  'w'  and  'P'  refer  to  values  at  the  wall  and  a 
point  P  close  to  the  wall.  The  point  yp+  is  established  by  the 
value  of  yp. 


(4.3.6) 


In  this  study  yp+  =15  for  most  cases. 


(iv)  Dissipation  rate  of  Turbulent  energy 


I  1 

ep  -  c/k//  (<yp)  (4.3.9) 

The  above  relations  describe  wall  functions  for  the  horizontal 
wall/  the  obvious  modifications  are  made  for  the  vertical  walls. 
4.3.2  Two-Pass  Procedure 

Figure  (5)  illustrates  the  solution  domains  for  each 
model  during  the  second  calculation.  In  implementing  the 
low-Fteynold  number  procedure  from  the  reattachment  point,  the 
calculation  proceeds  to  the  wall  with  grid  lines  placed  inside  the 
linear  sublayer-  This  creates  an  upstream  boundary  for  the 
low-Reynolds  number  solution  in  the  sublayer  at  the  vertical  plane 
through  the  reattachiuent  point.  Boundary  conditions  along  this 
plane  for  momentum  and  energy  in  this  region  are:  u+«y+  and  T+=y+ 
(as  discussed  in  the  previous  section).  The  boundary  condition  for 
the  turbulent  variable  k  is  obtained  hy  integrating  its  near-wall 
transport  equation  and  e  is  obtained  by  empirical  formula  directly 
fran  k,  as  outlined  in  the  previous  section. 

The  boundary'  conditions  for  the  entire  flow  field  during 
the  second  pass  are  sunmarized  in  Table  (4.1), 
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TABLE  4.1 

SUMMARY  OF  INLET  AND  BOUNDARY  OCNDITICKS 
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CHAPTER  V 


PREDICTIONS  OF  TURBULENT  REdRCLRATOG  BEAT  TRANSFER 
5.1  General  Remarks 

The  turbulent  models  described  in  Chapter  III,  together  with 
the  computational  procedures  of  Chapter  IV,  are  examined  by 
application  to  the  two-dimensional  back-step  and  the  pipe  expansion 
geometries  where  experimental  data  exist.  In  addition  to  the  heat 
transfer  experiments  of  Aung  and  Goldstein  [8] ,  Zemanick  and 
Dougall  and  Sparrow  and  O'Brien  [41]  (which  did  not  include 
detailed  hydrodynamic  measurements) ,  flow  field  comparisons  are 
made  with  the  pure  hydrcdynanic  measurements  of  other  investigators 

r  «sni 

All  calculations  are  performed  with  air  as  the  working  fluid. 
The  values  of  the  fluid  properties  of  air,  such  as  density  (  P  ) , 
dynamic  viscosity  (  u  ) ,  specific  heat  at  constant  property  (Cp)  and 
molecular  Prandtl  rnxnber  (Pr)  are  taken  from  Eckert  and  Drake 
[115] .  In  the  back-step  experiments  of  Aung  and  Goldstein  [8] ,  the 
walls  are  heated  to  uniform  tenperatures.  Ihe  difference  between 
free  stream  and  wall  tenperature  is  18°c.  For  the  pipe  expansion 
of  Zemanick  and  Dougall  [39] ,  the  wall  temperature  distribution  is 
given  in  Zemanick  [116] .  For  the  pipe  expansion  results  of  Sparrow 
sind  O'Brien  [41]  obtained  from  napthalene  sublimation  measurements, 
the  wall  tenperature  is  taken  as  the  sublimation  tenperature  of 
napthalene.  For  the  back  step  calculation  the  Reynolds  number  is 
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based  on  the  mean  inlet  velocity  and  step  height  and  for  the  pipe 
expansion  it  is  based  on  mean  inlet  velocity  ana  downstream  pipe 
diameter. 

Figure  (4)  shews  the  schematics  of  the  two-dimensional 
back-step  and  sudden  pipe  expansion  geometries  considered  in  this 
study.  The  point  y  »  0,  for  the  back-step  is  at  a  location  well 
outside  the  approach  boundary  layer,  and  for  the  pipe  expansion  it 
is  located  at  the  symmetry  axis.  Table  (5.1)  lists  the  actual 
dimensions  used  in  the  various  numerical  and  physical  experiments. 
In  this  table  all  distances  are  in  centimeters,  the  distance  x/H 
(in  tentis  of  step  or  expansion  height)  represents  the  final  station 
at  which  measurements  were  taken  and  the  expansion  ratio  d/D,  for 
the  back-step,  represent  the  actual  dimensions  of  the  upper  grid 
boundary  used  in  the  present  lumericai  study. 

The  results  presented  in  this  chapter  are  obtained  by  solving 
the  system  of  equations  of  Chapter  III  by  the  numerical  procedures 
of  Chapter  IV.  For  a  better  appreciation  of  the  tasks  involved  in 
the  present  study,  initial  calculations  are  presented  (mainly  for 
the  back-step  geometry)  just'fying  the  need  for  the  various 
improvements  and  modifications  implemented  during  the  development 
of  this  research.  Following  this  overview,  detailed  flow, 
temperature  and  turbulence  intensity  fields  and  heat  transfer 
calculations  obtained  with  the  two-pass  procedure  incorporating 
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T&BLE  5.1 

PHYSICAL  DIMENSIONS  CONSIDERED  IN  TEST  CASES 
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improvements  to  the  nodel  and  code  (as  discussed  earlier)  are 
presented.  Table  (5.1)  lists  the  geometries  considered  for  which 
comparisons  between  the  present  numerical  calculations  and 
experimental  measurements  are  made. 

Before  proceeding  directly  to  the  presentation  of  the 
numerical  results,  a  few  general  comments  regarding  the  calculation 
procedures  are  presented.  For  most  of  the  calculations  in  the 
standard  k-e  model  y£  =  15.  For  a  single  (first)  pass  with 
quadratic  differencing,  a  typical  non-uniform  grid  distribution 
used  is  (25  x  28),  convergence  is  obtained  after  about  250 
iterations  with  central  processor  times  on  an  IBM  3033  Processor 
Complex  averaging  about  4  minutes.  In  the  second  pass,  in  the 
redeveloping  region  the  number  of  horizontal  grids  are  increased  to 
about  32.  In  early  calculations,  using  hybrid  differencing  with 
the  standard  k-e  model,  a  typical  non-uniform  grid  used  was  (28  x 
32) ,  and  convergence  was  obtained  after  about  190  iterations.  Grid 
independence  of  results,  in  the  present  study  was  verified  for  all 
test  cases.  Figure  (9)  illustrates  the  verification  procedure  for 
the  back-step  experiments  of  Ajng  and  Goldstein  [8] .  Nusselt 
number  at  one  location  is  plotted  as  a  function  of  y^-grid 
densities.  The  x/H  location  chosen  is  purely  arbitrary,  for  most 
cases  a  value  just  downstream  of  the  reattachment  point  is  chosen. 
Hie  first  point  on  Figure  (9)  corresponds  to  15  y-grids  uniformly 
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spaced  across  the  channel.  All  the  other  points  correspond  to 

non-unlforra  grid  densities,  with  a  higher  concentration  of  grids  in 

the  near-wall  region.  As  the  number  of  horizontal  grid  densities 

are  increased  from  20  y-grids  to  about  35  y-grids ,  but  keeping  the 

near-wall  grid  location  the  same  so  as  to  maintain  the  v* 

V 

condition,  the  Nusselt  number  increases  slightly.  Beyond  30 
y-grids  the  changes  is  less  than  about  2  percent.  In  most  of  the 
present  calculations  the  number  of  y-grids  used  is  about  30.  A 
similar  verification  procedure  was  also  adopted  tor  the  x-grid 
densities. 

'Typical  grid  distributions  used  in  the  present  study  for  both 
the  back-step  and  pipe  expansion  geometries  are  shewn  in  Figures 
(7)  and  (10),  for  a  single  and  double  pass,  respectively.  The  grid 
distribution  shown  in  Figure  (7)  is  arranged  such  that  for  the 
entire  calculation  domain  the  near-wall  nodes  are  outside  the 
viscous  sublayer.  Typical  near-wall  distances  are  shown  in  the 
figure.  This  procedure  is  necessary  (as  discussed  previously)  when 
using  the  standard  k-e  model  in  a  single  pass.  Upon  establishing 
the  reattachment  length  by  the  previous  single-pass  procedure,  a 
typical  grid  distribution,  as  that  shown  in  Figure  (10),  is  then 
enployed  for  the  second  pass.  The  figure  also  shews  near-wall 
nodes  outside  the  viscous  sublayer  for  the  side  wall  and  up  to  the 
reattaciraent  point,  xRr  for  the  bottan  wall,  again  consistent  with 
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the  computational  domain  for  the  standard  k-e  model.  Beyond 
reattachment,  at  the  very  next  downstream  node,  grid  points  are 
placed  close  to  the  wall,  inside  the  viscous  sublayer,  and  the 
solution  procedure  utilizes  the  Low-Reynolds  number  version  of  the 
k-e  model.  The  initial  and  boundary  conditions  described  in  Table 
(4.1)  are  applied  at  the  following  locations:  initial  conditions  at 
the  plane  x  *  0,  boundary  conditions  at  the  near  wall  nodes,  free 
stream  (or  symmetry)  conditions  at  the  plane  y  =  0,  exit  conditions 
along  the  plane  x  =  NI  and  upstream  condition  for  the  Low-Reynolds 
number  model  along  the  plane  x  =  x^. 

A  brief  discussion  concerning  the  dependency  of  reattachraent 
length  with  Reynolds  number  and  step  height  is  now  presented  for 
the  back-step  geometry.  Fbr  turbulent  flows  reattachment  length  is 
independent  of  Reynolds  number.  Bowever,  for  laminar  and 
transitional  flews  there  is  a  significant  Reynolds  number 
dependency.  Initially  (laminar  regime),  reattachment  length 
increases  sharply  with  an  increase  in  Reynolds  number,  and  reaching 
a  maximum  (in  the  transitional  regime)  then  gradually  decreasing  to 
a  fairly  constant  value.  Figure  (11)  illustrates  the  effects  of 
Reynolds  number  (based  on  manentura  thickness  at  the  step),  on 
reattachment  length  for  laminar,  transitional  and  turbulent 
regimes.  Hie  experimental  data  points  are  that  of  Eaton  and 
Johnston  120 J  and  the  solid  line  joins  the  experimental  points. 
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Initial  calculations  were  performed  using  three  different  step 
heights  to  examine  the  effects  of  step  heights  on  reattachment 
lengths.  Figure  (12)  illustrates  that  Xp/H  remains  constant  with 
H.  These  results  were  obtained  by  using  the  original  standard  k-e 
model  (i.e.  no  modifications  to  <^),  However,  the  same  constant 
behavior  is  expected  when  using  the  inproved  model,  the  only 
difference  lies  in  the  magnitude  of  the  reattachment  length. 
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5.2  Initial  Results  (Single  Pass) 

Initial  calculations  were  performed  for  the  two-dimensional 
back-step  and  pipe  expansions  by  incorporating  both  the  mixing 
length  and  standard  k-e  models.  For  the  k-e  model  was  assumed  a 
constant  value  of  0.09  and  Prt  approximated  fcy  Bqn.  (3.4.2). 
Approach  and  boundary  conditions  used  have  been  discussed 
previously.  Presented  here  is  a  summary  of  typical  initial  heat 
transfer  results  for  the  back-step  in  which  Nusselt  number 
(normalized  by  its  value  just  upstream  of  the  step)  is  plotted  at 
various  downstream  distance  (normalized  by  step  height)  for  a 
Reynolds  number  (based  on  step  height)  of  1728.4.  More  extensive 
standard  k-e  results  obtained  in  the  present  research  including 

4m 

hydrodynamic  and  turbulence  flow  field  predictions  are  presented  in 
Ref.  1106 J .  In  results  to  be  presented  here  for  the  initial 
calculations  comparison  is  made  between  numerical  prediction  and 
experimental  data  of  Aung  and  Goldstein  [8J  for  their  highest 
Reynolds  number  case.  Figure  (13)  shows  results  obtained  frari  the 
initial  calculations  in  which  the  mixing  length  model  was  employed. 
It  is  apparent  from  the  figure  that  this  sirple  model  consisting  of 
a  single  turbulence  scale  does  not  account  for  the  canplex 
turbulence  phenomena,  particularly  in  the  recirculating  region. 
This  is  evident  by  the  comparitively  small  values  of  heat  transfer 
obtained  when  using  the  mixing  length  model.  Figure  (14)  shows 
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heat  transfer  results  obtained  for  the  same  case  by  the  standard 
k-e  model.  With  this  model,  the  predictions  consistently 

overestimate  the  experimental  data  particularly  in  the  redeveloping 
region.  However,  the  overall  agreement  is  fairly  satisfactory 
considering  the  inherent  ’simplicity’  of  the  k-e  model.  The 

constraints  on  the  choice  of  the  near  wall  grid  line  in  these 

calculations,  as  discussed  previously,  are  that  it  must  lie  outside 
the  linear  sublayer  and  yet,  be  close  enough  to  the  wall, 

particularly  for  heat  transfer  calculations ,  to  provide  accurate 
resolution  of  the  flow  in  the  near^-wall  region.  For  most  of  the 

present  calculations  y4"  averages  about  15,  which  is  the  approximate 

P 

average  value  in  the  recirculation  zone. 

Since  the  Reynolds  nimbers  of  the  present  flow  calculations 
for  the  step  are  comparatively  low,  the  physical  distance  of  y£ 
from  the  wall  is  large  and  therefore  resolution  of  the  details  of 
the  near-wall  flew  is  not  possible.  In  an  attempt  to  improve  the 
results  by  increasing  resolution  near  the  walls  the  IcvHReynolds 
number  form  of  the  model  was  employed  for  the  entire  region 
downstream  of  the  step.  A.  typical  result  obtained  by  utilizing 
this  model  is  shown  in  Figure  (15).  The  heat  transfer  results 
underestimate  the  data  for  the  recirculating  region,  but  fairly 
good  agreement  is  obtained  for  the  redeveloping  region.  This 
behavior  is  understandable  since  the  justifications  of  the 


low-Reynolds  number  terms  in  this  model  was  based  on  studies  with 
accelerated  boundary  layer  type  flows  L 60] .  Such  a  model  was 
needed  for  it  yielded  an  additional  improvement  by  permitting 
resolution  of  the  near-wall  region  downstream  of  reattachment. 
This  is  extremely  important  when  a  local  extremum  in  the  velocity 
profile  occurs  near  the  wall  as  is  the  tendency  in  the  present 
cases.  The  results  obtained  by  the  standard  k-e  model  and  the 
low-Reynolds  number  model  indicate  that  the  standard  model  gives 
better  overall  prediction  in  the  recirculation  zone  and  that  the 
low-Reynolds  number  model  is  better  for  the  redeveloping  region. 
The  two-pass  procedure  exploits  this  phenomena.  However.-  the 
feasibility  of  such  a  procedure  would  depend  on  the  reliability  of 
predicting  reattachment  lengths.  This  is  accomplished  by 
incorporating  various  improvements  in  the  standard  k-e  model 
especially.  Rssults  utilizing  these  improved  models  are  discussed 
in  the  next  section. 

To  yield  overall  better  predictions  of  heat  transfer  results 
another  variation  of  the  model  was  tested.  This  is  discussed  in 
Section  (3.4.3)  with  the  additional  term  to  the  e -equation  given  by 
Bqn.  (3.4.9) .  This  procedure  is  often  referred  to  as  the 
'sensitivity'  correction  of  Hanjalic  and  Launder  162] .  Results 
obtained  by  this  procedure  are  shewn  in  Figure  (16).  Again,  no 
substantial  improvement  in  heat  transfer  results  are  obtained. 
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Predictions  of  flow  field  variables  by  the  two-pass  procedure 
will  be  presented  in  a  later  section.  All  of  the  above 
calculations  (except  for  the  mixing  length  model)  yield 
reattachment  lengths  of  about  4.5  step-heights  downstream  of  the 
step,  with  the  maximum  (peak)  value  of  Nusselt  number  occuring 
about  one  step-height  upstream  of  reattachroent.  (Reattachment 
lengths  obtained  by  the  various  procedures  are  surrmarized  in  the 
next  section).  The  value  of  4.5  is  consistent  within  the  range  of 
values  obtained  experimentally  by  Aung  and  Goldstein  18 J  for  a 
tripped  approach  boundary  layer.  However,  the  models  in  their 
present  form  underestimate  the  reattachment  length  for  all  the 
other  test  cases  considered.  This  phenomenon  is  also  discussed  in 
detail  in  the  next  section. 


11U 


5.3  Improved  Results  (Single  Pass) 


As  discussed  in  Qiapter  III  the  standard  k-e  model  and  the 


low-Reynolds  number  form  of  the  model  were  iirproved  in  the  present 
work  by  functionalizing  and  Prt  through  Bgns.  (3.4.36)  and 
(3.4.37),  respectively.  These  (and  all  subsequent)  calculations 
incorporated  quadratic  approximations  for  the  convective  terms,  as 
discussed  previously.  Figure  (17)  conpares  the  heat  transfer 
coefficients  obtained  by  the  quadratic  and  by  hybrid  differencing 
schemes.  The  quadratic  procedure  yields  slightly  iirproved 
predictions  of  heat  transfer  for  the  recirculating  region. 
However,  improvements  in  predictions  of  reattachment  lengths  are 
more  significant  and  further  the  scheme  requires  lower  near-wall 
grid  densities  for  grid  independent  solutions,  as  indicated  on  the 


figure. 

Figure  (18)  conpares  heat  transfer  obtained  by  both  the 
present  unproved  standard  k-e  and  low-Reynolds  number  models.  The 


trend  is  the  s«one  as  discussed  previously.  For  instance,  the 
standard  model  and  low-Reynolds  number  form  of  the  model  better 
represent  the  flow  in  the  recirculating  and  redeveloping  regions, 
respectively,  thus,  re-eitphasising  the  need  for  the  two-pass 
procedure.  Figure  (19)  shows  similar  heat  transfer  results,  but 


for  a  different  value  of  Reynolds  number. 

In  an  attempt  to  yield  better  overall  heat  transfer 
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predictions  for  the  entire  region  by  the  improved  standard  k-e 
model/  without  resorting  to  the  two-pass  procedure  for  which 
results  are  described  later,  the  ’sensitivity’  correction  of 
Leschziner  and  Rodi  164J  were  tested.  In  this  procedure  to  the 
e -equation  an  additional  term  is  introduced,  as  given  by  Bgn. 
(3.4.10).  Results  obtained  by  this  procedure  are  shorn  in  Figure 
(20).  Again,  no  substantial  improvement  in  heat  transfer  results 
for  the  redeveloping  region  are  obtained. 

Heat  transfer  results  for  the  pipe  expansion  case  (for  a 
Reynolds  number  of  66,260)  fcy  the  improved  standard  k-e  model  is 
presented  in  Figure  (21).  Also  indicated  on  the  figure  are  the 
experimental  data  of  Zemanick  and  Dougall  L39] .  Again  predicted 
results  are  larger  than  experimental  ones.. 

The  two-pass  procedure  described  in  Chapter  III  requires 
partitioning  of  the  flow  field  into  two  regions:  recirculating  and 
redeveloping  regions.  Consequently,  the  applicability  of  the 
procedure  relies  on  the  accuracy  and  consistency  to  which  the 
reattachment  length  can  be  predicted.  It  is  for  this  reason  that 
various  numerical  experiments  were  undertaken  with  emphasis  on 
reliable  predictions  of  reattachment  lengths.  Table  (5.2) 
summarized  reattaclment  lengths  obtained  by  the  various  turbulent 
approximation  procedures  discussed  previously.  In  the  table,  the 
expression  "partially  improved  C^"  refers  to  exclusion  of 
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wall-danping  corrections  and  imposing  the  equilibrum  of  turbulence 
energy  assumption  in  the  derivation  of  the  function  for  Cy,  i.e. , 
A"  «  0  and  P/e  =  1  in  Bgn.  (3.4.36).  Fully  improved  C^,  refers  to 
Egn.  (3.4.36).  The  results  indicate  that  by  employing  the  fully 
inproved  expression  for  Cy  in  the  standard  k-e  model  (together  with 
quadratic  differencing  for  the  convective  terms)  yields 
reattachment  lengths  that  are  consistent  with  experiments,  except 
for  the  Aung  and  Goldstein  18]  case.  However,  the  reason  for  the 
over-prediction  is  explainable.  In  their  experiments  the  boundary 
layer  trip  was  located  at  eleven  trip-heights  upstream  of  the  step, 

i  nK  vopmI  f  e<  ■»  h  a  Ki’mV»q>*  1  atu^I  /vF  furhtil  oniNa  i  nfapc  i  Kf  af  fV»« 
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step.  In  order  to  gain  same  appreciation  for  the  magnitude  of  this 
effect,  a  numerical  experiment  was  conducted  by  using  the  inproved 
k-e  model  to  perform  a  turbulent  back-step  flew  field  calculation 
with  the  step  height  taken  as  the  height  of  the  trip.  Tne  results 
then  were  examined  to  obtain  quantitative  estimates  of  the  level  of 
k  and  e  at  11  trip  heights  downstream.  These  values  were  then  used 
as  upstream  condition  input  for  the  regular  back-step  calculation 
in  lieu  of  the  method  described  earlier. 

Although  the  effect  on  heat  transfer  of  the  higher  level  of 
turbulence  is  slight,  as  shewn  in  Figure  (22),  the  reattachraent 
length  is  reduced,  and  essentially  ccnparable  with  the  experiment. 
The  reattachment  length  obtained  is  indicated  in  the  last  column  of 
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Table  (5.2).  Since  there  remains  sane  questions  of  whether  or  not 
the  physics  are  really  iirproved  by  this  method  of  obtaining 
upstream  conditions  and  also  as  to  whether  or  not  heat  transfer 
results  using  it  are  actually  in  closer  agreement  with  experiment, 
it  was  not  adapted  for  the  final  calculations. 
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5.4  Improved  Results  (Two-Pass) 

The  two-pass  method  was  the  procedure  which  as  previously 
indicated  gave  the  best  overall  prediction  of  heat  transfer.  In 
applying  it  during  the  first  pass,  calculations  are  performed 
according  to  the  previous  section  to  obtain  the  reattachment 
length.  Specifically,  for  the  best  results,  the  first  pass 
employed  the  improved  standard  k-£  model,  with  and  Prp  given  by 
Bqns.  (3.4.16)  and  (3.4.37),  respectively,  fbr  the  second  pass, 
the  above  procedure  is  maintained  up  to  reattachnvent,  and  beyond 
reattachnvent  the  improved  lcw-Reynolds  number  form  of  the  model  is 
used  with  C^,  and  Prt  given  by  Bgns.  (3.4.38)  and  (3,4*39), 
respectively. 

5.4.1  Results  for  Back-Step 

Before  presenting  the  heat  transfer  predictions  frcnv  the 
two  pass  method,  the  hydrodynamic  and  turbulence  flow  field 
predictions  are  discussed.  The  purpose  is  to  illustrate  that 
although  the  enphasis  of  the  present  study  is  on  heat  transfer,  the 
present  procedure  reliably  predicts  the  other  turbulent  flew  field 
variables  as  well.  Typical  computed  velocity  profiles  for  two 
different  experimental  cases  are  shown  in  Figure  (23)  and  (24).  As 
can  be  seen,  the  agreement  between  the  experimental  data  and 
present  predictions  is  good.  Figures  (25)  and  (26)  represent 
cross-stream  pressure  variation  and  surface  pressure  coefficient 
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downstream  of  the  step,  respectively.  In  the  recirculating  region, 

the  pressure  is  not  uniform,  as  is  expected,  but  approaches  a 

uniform  value  in  the  redeveloping  boundary  layer.  Figure  (27)  is  a 

comparison  of  ccnputed  temperature  profiles  for  one  Reynolds  number 

with  the  corresponding  data  of  Aung  and  Goldstein  (8J  and  is 

typical  of  the  agreenent  of  results  at  the  other  two  Reynolds 

number  as  well.  Figure  (28)  and  (29)  compare  mean  velocity  and 

temperature  profiles  in  the  reattached  boundary  layer  with  their 

universal  'log-wall'  functions.  ihe  ccnputed  slopes  deviate 

slightly  fron  that  of  the  universal  slopes.  This  is  typical  of  a 

reattached  boundary  layer  18,25].  Figures  (30)  and  (31)  depict  the 

ccnputed  behavior  of  turbulent  kinetic  energy.  Figure  (30)  shows 

the  profiles  at  selected  streamwise  locations;  the  experimental 

2 

data  points  represent  u'  only.  Hie  same  qualitative  trend  is 

% 

observed  in  the  figure  between  the  prediction's  and  the  experiment, 

and  further  in  the  recirculating  region  k  is  about  40  percent 

,2 

higher  than  u  [25] ,  since  in  this  region  all  three  components  of 
the  Reynolds  stresses  are  inportarit.  Figure  (31)  is  a  plot  of  the 
streamwise  variation  of  the  maximum  (peak)  value  of  k  and  the 
experimental  data  points  are  the  total  k  values.  As  can  be  seen 
that  the  agreement  between  predictions  and  experimental  data  is 
again  very  good. 

A  comparison  of  the  streamwise  variation  of  the  ccnputed 
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heat  transfer  coefficients  with  the  experimental  data  of  Aung  and 
Goldstein  [33  is  given  in  Figures  (32)  through  (34).  The  overall 
agreement  is  acceptable.  Figure  (32)  sirmarizes  the  results 
obtained  by  the  improved  procedure  with  a  single  pass  only  and  the 
two-pass  procedure. 

5.4.2  Results  for  Pipe  Expansion 

Heat  transfer  results  for  three  pipe  expansion  cases  are 
plotted  in  Figures  (35)  through  (37).  The  results  of  previous  k- 
calculations  by  Chieng  and  Launder  159  J  together  with  the 
experimental  data  of  Zemanick  and  Dougall  L39J  are  also  shewn  for 
comparison.  In  the  recirculating  region  the  llusselt  number 
obtained  by  the  present  improved  standard  k-e  model  differs 
significantly  from  that  obtained  by  Chieng  and  Launder  (59J .  The 
reasons  for  this  difference  are  that  Chieng  and  Launder  15SJ  in 
their  calculations  assumed  constants  for  both  Cy  and  Pr^  and  more 
importantly,  their  calculations  proceeded  from  the  buffer  zone. 
Starting  the  solution  inside  the  buffer  zone  is  only  possible  by 
specifying  k  and  e  for  this  region.  Early  attempts  were  made  in 
the  present  study  to  incorporate  the  Chieng  and  launder  (59J 
functions  valid  for  the  buffer  zone  11U6J .  However,  for  sane 
unresolved  reason  the  predictions  obtained  did  not  correspond  with 
th.:  results  of  Chieng  and  Launder  L59J .  It  is  obvious  t>iat  the 
results  from  the  present  two-pass  procedure  are  a  narked 
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improvfsient  over  the  Chieng  and  Launder  [59]  calculations  and  also 
are  in  reasonably  good  agreement  with  experimental  results  in  the 
entire  flow  field. 

Figure  (38)  shows  the  maximum  (peak)  Nusselt  number  dependence 
on  Reynolds  number.  The  predictions  exhibit  a'  slightly  greater 
sensitivity  to  Reynolds  number  than  do  the  measurements.  This 
slight  difference  is  insignificant  when  considering  the  range  of 
Reynolds  number  studied.  The  predicted  maximum  Nusselt  number  is 
represented  by  the  following  expression. 

=  c*eH°‘75'  c  “  °*145 
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5.5  Downstream  Face  Heat  Transfer 

Figure  (39)  shows  the  side-wall  heat  transfer  coefficient  for 
the  back-step  for  two  Reynolds  numbers.  The  Nusselt  nutter  is 
normalized  by  its  average  value  and  the  horizontal  axis  is 
normalized  to  vary  from  U  to  1,  the  former  corresponding  to  a  point 
adjacent  to  the  opening  and  the  latter  to  a  point  where  the 
side-wall  meets  the  horizontal  wadi,  ihe  results  obtained  (tiicugh 
not  possibly  evident  from  the  figure)  indicate  that  the  Nusselt 
numbers  are  very  small  compared  to  say  the  base  wall  (about  10 
percent  of  base-wall  value).  The  highest  heat  transfer  occurs  on 
the  portion  of  the  side-wall  adjacent  to  the  main  flaw  and  the 
lowest  occurs  in  the  comer  region.  No  experimental  data  are 
available  for  ccnparison. 

Figure  (40)  through  (42)  shew  similar  predicted  results  of 
heat  transfer  coefficients  for  three  Reynolds  nunber  for  the 
downstream  face  of  a  pipe  expansion.  Carparison  is  made  with  the 
experimental  data  of  Sparrow  and  O'Brien  [41].  Agreement  with  the 
experiment  is  better  for  the  larger  Reynolds  numbers.  This 
observed  phenomena  can  be  explained  by  the  fact  that  in  the  region 
adjacent  to  the  downstream  facing  wall  the  improved  standard  k-e 
model  is  used;  and,  as  mentioned  previously,  the  near-wall 
resolution  of  this  model,  particularly  for  low-Reynolds  nunbers,  is 
severly  restricted.  Figure  (43)  shews  the  dependence  of  the 
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a  veraye  Nusselt  number  on  Reynolds  number,  the  growth  law  is 
essentially  the  same  as  that  of  the  experimental  data.  The 
predicted  growth  law  is  represented  by  the  following  expression. 
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5.b  Correlations  of  Back-Step  Heat  Transfer 

Results  related  to  correlations  of  heat  transfer  coefficients 
are  presented  in  Figures  (44)  through  (48)  for  the  bade  step 
geometry*  All  calculations  vrere  performed  by  assuming,  at  the 
step,  an  untripped  turbulent  boundary  layer  having  a  heated 
starting  length.  Also,  the  same  momentum  thickness  for  the 
incident  boundary  layer  is  assigned  to  all  Reynolds  numbers.  This 
is  not  too  damaging  an  assumption  because  the  heat  transfer 
coefficients  beyond  the  step  are  fairly  independent  of  approach 
boundary  layer  thicknesses.  This  is  demonstrated  in  Figure  (44) 
for  three  different  approach  boundary  layer  thicknesses.  The 
region  that  shews  the  greatest  sensitivity  of  heat  transfer 
coefficients  to  different  approach  boundary  layer  thicknesses  is 
the  region  ^ordering  the  reattachroent  point. 

Figure  (45)  is  a  logarithmic  plot  of  Nusselt  number,  Nuj, 
based  on  distance  x,  downstream  of  reattachraent  versus  Reynolds 
number,  Rej#  for  five  Reynolds  nunber  cases.  The  flat  plate 
turbulent  boundary  layer  theory  for  heat  transfer  is  also  shown  for 
comparison.  The  computed  results  all  tend  to  the  0.8  power  growth 
law  of  the  flat  plate  relation  sufficiently  far  downstream  of  the 
step.  Here  this  occurs  between  10  to  13  step  heights  downstream  of 
the  step.  These  results  demonstrate  the  fact  that  the  flat  plate 
boundary  layer  theory  based  on  distance  downstream  of  reattachment 
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cannot  be  used  to  predict  the  heat  transfer  in  the  redeveloping 
boundary  layer  region.  Even  at  significant  distances  downstream 
where  the  sane  power  law  behavior  is  followed,  it  will  always  tend 
to  underestimate  the  heat  transfer  rate.  In  the  present 
predictions,  the  confuted  results  exceed  the  flat  plate  theory  by 
factors  ranging  from  1.3  to  1.8. 

Figures  (46)  and  (47)  show  the  dependence  of  the  maximum  (peak) 
and  side-wall  Nusselt  numbers  with  Reynolds  numbers,  respectively. 
This  power  law  dependence,  for  convenient  reference  purposes,  is 
summarized  in  Table  (5.3)  and  should  be  of  use  to  the  designer. 
Hie  power  law  behavior  for  the  maximum  Nusselt  number  is  used  in  an 
attenpt  to  correlate  the  distribution  of  heat  transfer  coefficients 
for  design  purposes.  This  result  is  presented  in  Figure  (48). 
Figure  (48)  presents  results  for  Nu  Re  0*S8  versus  x/H,  where  the 
exponent  0.68  represents  the  power  law  dependence  of  the  maximum 
Reynolds  number.  Also  plotted  on  the  figure  are  the  experimental 
data  of  Aung  ana  Goldstein  [8] .  The  data  points  fall  within  the 
expected  range,  except  for  the  peak  values.  This  is  probably  due 
to  the  slight  influence  of  the  difference  in  approach  boundary 
layer  thicknesses  on  heat  transfer  coefficients.  The  graphs  of 
Figures  (46)  through  (48)  can  be  used  with  some  confidence  to 
predict  back-step  heat  transfer  coefficients  with  accuracy 
sufficient  for  most  design  purposes,  within  the  range  of 
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TABLE  5,3 

CORRELATION  OF  MAXIMUM  (PEAK)  AND  AVERAGED  (FOR  DOWNSTREAM  FACE) 
NUSSELT  NUMBER  WITH  REXNCDLS  NUMBER  FOR  BACK  STEP 


Reynolds  number  considered.  As  a  final  corment,  the  results  shown 
in  Figures  (46)  through  (48)  are  very  restrictive  in  the  sense  that 
they  represent  heat  transfer  coefficients  for  a  constant 
step-height  (H)  of  0.64  cm.  and  a  constant  approach  boundary  layer 
thickness  of  1.5H.  Future  studies  to  be  undertaken  will  examine 
the  effects  of  step-height  on  heat  transfer  coefficient  beyond  the 
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5.7  Discussion  of  Results 

The  present  results  indicate  that  using  the  usual  two-equation 
models  of  turbulence  (standard  and  low-Reynolds  number  form  with 
a  constant)  yields  unacceptable  predictions  of  reattachment  lengths 
and  heat  transfer  coefficients  for  the  entire  flow  field  downstream 
of  a  sudden  change  in  geometry.  A  great  improvement  in  results  are 
obtained  by  using  the  improved  versions  of  the  two-equation  models 
in  the  two-pass  procedure.  Comparison  of  computed  results  with 
experimental  results,  and  with  calculations  of  previous 
investigators  indicate  that  at  moderate  Reynolds  number  the 
two-pass  method  with  the  improved  two  equation  models  allows 
accurate  calculation1  of  heat  transfer  over  the  entire  flew  field. 


CHAPTER  VI 


CONCLUDING  REMARKS 

<j.l  Achievements  of  Present  Analysis 

The  results  presented  in  Chapter  V  demonstrate  that  by 
employing  modified  versions  of  the  k-e  model  applied  in  a  sequence 
of  two  passes,  predictions  of  heat  transfer  coefficients  can  be 
obtained  fairly  accurately  for  recirculating  and  subsequent 
redeveloping  flows  beyond  abrupt  expansions  where  the  iteynolds 
numbers  range  from  low  to  moderate  values.  The  procedures  developed 
herein  impose  minim  al  restrictions  on  near-wall  grid  densities, 
are  numerically  efficient  and  employ  fairly  simple  turbulence 
modeling.  They  are  therefore,  particularly  appropriate  for 
engineering  applications. 

Important  turbulence  phenomena,  such  as  streamline  curvature 
effects,  and  pressure-strain  (scalar)  interactions  including 
'wall-damping'  effects  have  been  incorporated  in  the  k-e  model 
through  functionalizing  two  parameters  (C^  and  Pr^)  that  were 
heretofore  assumed  as  constants.  Thus,  some  of  the  most  important 
turbulence  interactions  are  included  in  this  improved  k-e  model  and 
at  the  same  time,  side-stepping  the  cumbersome  numerical 
calculations  necessary  for  the  more  ccnplex  higher  order  schemes. 
This  functionalization  is  accomplished  by  employing  the  basic  ideas 
of  the  'second-order*  algebraic  stress  level  of  closure  without 
actually  employing  the  full  method  in  the  calculation. 
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Hie  eventual  modeling  and  numerical  procedures  adopted  for  the 
present  study  yield  consistently  acceptable  predictions  of 
reattachment  lengths,  correlations  of  heat  transfer  coefficients 
and  flow  and  turbulence  intensity  fields  of  the  entire  flow  field. 
Throughout  the  course  of  this  study,  enphasls  has  focused  cm 
evaluating  the  capabilities  and  feasibilities  of  both  the  numerical 
and  turbulent  closure  procedures.  Consequently,  the  numerical  code 
used  is  structured  in  such  a  way  that  it  can  be  readily  used  to 
evaluate  various  modifications  of  turbulent  viscosity  models. 
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6.2  Limitations  of  Present  Analysis 

'Hie  major  shortcoming  of  the  present  analysis  lies  in  the 
inherent  assunpr.ion  that  underlies  all  turbulent  viscosity  models, 
i.e.  use  of  the  Boussinesq  assumption.  By  this  assumption  the 
shear  stress  is  forced  to  instantaneously  respond  to  changes  in  the 
mean  strain  field.  This  concept  is  not  strictly  valid  in 
non-equilibruim  flews,  of  which  recirculating  flow  unfortunately  is 
an  example-  A  shortcoming  of  the  two-pass  procedure  used  is  that 
it  requires,  a  priori,  fairly  accurate  prediction  of  reattachnent 
lengths.  For  abrupt  expansion  gec*netries,  as  those  examined  in 
this  study,  fairly  accurate  reattachment  lengths  were  obtained. 
However,  for  other  engineering  heat  exchanger  geometries  it  is  not 
certain  if  the  same  level  of  success  can  be  achieved.  Uiere  are 
other  minor  limitations  of  the  present  research  which  become  even 
less  significant  when  viewed  in  the  context  of  applications  of  heat 
exchange  equipment;  for  example,  neglect  ion  of  the  intermittent  and 
fine  scales  structure  of  turbulence.  As  a  final  renark,  Chandrsuda 
and  Bradshaw  L25]  presented  detailed  measurements  of  Reynolds  shear 
stress  and  triple  products  for  the  entire  flow  field  beyond  the 
back-step.  Iheir  results  indicate  that  the  triple  products  play  an 
important  role  in  the  ’spreading’  of  the  shear-layer,  particularly 
around  reattachnent.  However,  the  k-e  viscosity  model  in  its 
present  form  is  not  ’sensitized’  to  reflect  triple  products. 
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6.3  Suggestions  for  Future  Studies 

(a)  Direct  application  of  the  procedures  developed  in  this 
study  can  be  applied  to  study  turbulent  heat  transfer 
flows  over  a  cavity,  for  which  experimental  measurements 
are  reported  146J .  For  such  a  geometry  the  flow 
recirculates  inside  the  cavity  and  the  original  shear 
layer  is  only  slightly  disturbed,  consequently,  the 
single  pass  procedure  might  suffice. 

(b)  By  making  the  appropriate  changes  in  the  transport 

equations,  buoyancy  driven  flews  (such  as  that  in  a 
vertical  sudden  expansion)  can  be  studied.  It  is 

expected  that  the  present  improvements  to  the  code  and  to 
the  basic  model  are  adequate  for  analyzing  such  flow 
fields. 

*•'  «» 

(c)  Further  improvements  on  modeling  the  dissipation  equation 
is  necessary.  The  inproved  model  should  more  precisely 
account  for  near-wall  effects  and  probably  triple 
products.  If  such  a  model  is  devised,  then  it  is 
possible  to  predict  results  for  the  entire  flow  field  by 
a  single  pass. 

(d)  In  present  studies  the  length  scales  (kV^/e)  for  both 
the  velocity  arid  tenperature  are  assured  proportional. 
As  discussed  previously  this  will  only  be  true  if  P/e  « 
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Pg/Ep,  which  is  rarely  approached  in  a  recirculating 
region.  The  way  forward  would  be  to  replace  the 
temperature  length  scale  by  0*2 >3/2^  and  then  obtain 

differential  equations  of  transport  for  both  e'2  and  e'. 

0 


APPENDIX  A 


CLOSURE  LEVELS  OF  TURBULENCE 

lhe  mathematics  of  same  of  the  most  canton  closure  techniques 


for  the  stresses  u^u!  and  uje'  are  presented  in  this  appendix. 

(1)  Turbulent  Viscosity  Models 

According  to  the  Boussinesq  [54]  idea,  the  tensorial 
representation  of  the  Reynolds  stresses  and  scalar  fluxes, 
respectively,  are. 
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where  ut  and  Ft  are  the  turbulent  viscosity  and  diffusivity  of 

,  I  I 

temperature  respectively.  Tie  term  1/3  6..  u  u  in  Ban.  (Al) 

*  ij  mm 

simply  ensures  that  the  norm]  stresses  do  not  vanish  in  the 
absence  of  mean  strain,  and  they  are  equal  in  magnitude.  The 
approximations  for  the  turbulent  viscosity,  will  now  be 

discussed. 


(a)  Algebraic  Expressions 

By  Prandtl's  mixing  length  hypothesis  [53], 
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where  lQ  is  the  mixing  (characteristic)  lengtn,  and  is 
approximated  as. 


*0  =  KY 
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where  k  is  Von  Karman  constant  {=0.4).  Such  an  approximation  for 
£q  is  clearly  only  valid  for  regions  remote  from  the  wall, 
excluding  the  sublayer  and  buffer  layer  close  to  the  wall.  An 
improvement  of  this  length  scale  was  suggested  by  Van  Driest  156] 
for  the  near-wall  region.  The  expression  is 

*  k  y  |l  -  exp  (-y+  xw/  A+)  J  (AS) 

In  the  fully  turbulent  region  Eqn.  ( A5)  reduces  to  Bqn.  (A4) ,  as  is 
expected. 

(d)  Differential  fororessions 

In  this  procedure  the  viscosity  is  again  approximated  by 
scales  (length  or  velocity)  which  are  in  turn  approximated  by 
transport  differential  equations. 

(i)  One  Equation  Model 

This  modeling  procedure  has  been  developed  and  tested  by 
Bradshaw  et  al  [57] .  The  viscosity  is  approximated  by 

1  _ 

Pt  =  C  p  ,  k  =  ~  (uiui)  (A6) 

where  C  is  a  constant,  k^  is  the  velocity  scale  and  (as 
described  above)  is  the  length  scale.  The  turbulent  kinetic 
energy,  k,  is  expressed  by  a  transport  differential  equation.  The 
differential  equation  is  obtained  by  a  tensorial  contraction  of  the 
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the  Reynolds  stress  equation.  The  Reynolds  stress  equation  is 
developed  in  the  next  section.,  thus,  development  of  the 
differential  equation  for  k  is  deferred  to  the  next  section. 

(ii)  Two  Equation  Model 

This  modeling  procedure  is  bated  on  approximating  the 
turbulent  viscosity  by  two  scales,  length  and  velocity  [60 J . 


where 
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and  Cu  normally  assumed  a  constant  value  of  0.09.  The  dissipation 
of  turbulent  energy,  e  ,  is  expressed  by  a  transport  differential 
equation.  For  the  same  reason  as  mentioned  for  k,  the  e  equation 
is  deferred  to  the  next  section. 

(2 )  Stress  and  Flux  Transport.  Models 

Transport  equations  for  the  Reynolds  stress  and  scalar 
flux  is  derived.,  and  a  brief  discussion  of  the  closure  tectintques  is 
discussed,  (The  unsteady  form  of  the  governing  equations  are 
saintained,  although  in  the  present  study,  the  steady  tom  of  the 
equations  are  considered).  The  instantaneous  values  of  \he 
Navier-Stokes  equation,  Bqn.  (1.2.1)  are  replaced  by  the  sun  of  the 


mean  and  fluctuating  quantities,  as  described  in  Section  {1.2),  and 


applying  the  time-averaged  continuity  equation  to  further  sinplify, 
results  in  a  differential  equation  for  the  fluctuating  velocity. 
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Writing  a  similar  equation  for  the  uj  fluctuation,  and  multiplying 
the  uj  equation  by  uj,  and  the  uj  equation  by  uj,  and  adding  the 
two  resulting  equations,  yields, 
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Time-averaging  Bqn.  (AS)  and  grouping  sane  of  the  expressions 
results  in  the  Reynolds  stress  transport  equation  for  the 
correlation  uju7  Omitting  the  usual  overbars  for  the  mean  fluid 
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dynamic  variables  results  in  Bqn.  (3.4.13). 

A  similar  analysis  is  performed  for  the  scalar  (energy) 
equation,  Hie  procedure  is  as  follows:  substitute  the  Reynolds 
procedure  of  Section  (1.2)  into  the  instantaneous  energy  equation, 
Bqn.  (1.2.4),  and  multiply  the  new  equation  by  u^,  the  result  is 
then  added  to  the  product  of  Eqn.  (A8)  and  6'.  Time-averaging  the 
resulting  equation  yields  the  transport  equation  for  the  scalar 
flux,  u!  e'.  Again  by  emitting  the  usual  overbars  for  the  mean 
quantities,  Bgn.  (3.4.14)  is  obtained. 

The  transport  differential  equations  for  k  and  e,  as  used  in 
the  two-equation  model,  are  now  developed.  Performing  a 
contraction  on  the  Reynolds  stress  equation,  Bqn,  (3,4,13),  i.e., 
setting  i=j,  and  replacing  u!u ^/2  by  k  results  in  the  transport 
equation  for  k, 
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Eqn.  (AlO)  is  the  differential  equation  that  governs  the  transport 
of  turbulent  kinetic  energy,  k,  and  Terms  1,11,  and  IV  have  the 
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same  physical  interpretations  as  these  of  Egn.  (3.4.13). 

The  derivation  of  the  transport  equation  for  e  ,  the 
dissipation  rate  of  turbulent  energy,  is  slightly  more 
algebraically  involved.  The  mathematical  procedure  is  briefly 
outlined  below.  Ey  definition,  frcm  Bqn.  (3.4.13)  the  dissipation 
(Term  II) ,  e  ,  is,  ' 
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Differentiate  Bqn.  (A8;  with  respect  to  3/3  x£  and  multiply 
throughout  by  2  v3  Ui/3x  ,  time-average  the  results  and  use  the 
definition  of  e  ,  Bqn  (Ail).  By  grouping  various  terms,  the  result 
is  the  transport  equation  for  e  , 
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Ter™?  I, II,  and  III  retain  their  respective  physical  meanings  as 
given  in  Bqn.  (3.4.13).  Ihese  differential  equations  for  the 
transport  of  k  and  c  mist  be  closed  so  as  to  be  applicable.  The 


closure  of  Bgn.  (AlQ)  for  the  turbulent  kinetic  energy  is 
accomplished  by  comparing  each  of  the  terms  with  experimental 
measurements,  and  deriving  empirical  relationships.  Ft>r  the  e 
equation,  Bqn  (A12) ,  this  is  not  possible  for  none  of  the  terms  are 
accessible  to  measurements;  however,  closure  is  accomplished  by 
gross  simplifications.  Launder  169!  presented  the  closure  details 
for  both  equations. 

Within  the  framework  of  the  turbulent  viscosity  model,  Bqns. 
(All)  and  (Ai2)  reduces  to  [69J , 


where  o,  ,  o  ,  Qfrl  and  C  c  2  ars  empirical  constants,  their  values 

iC  G 

are  listed  in  Table  (3.1).  _ _ 

,2 

Often  analogous  transport  equations  are  needed  for  6  and 
eQ(«2  r  Oe'/3xk)2).  However,  in  this  study,  as  will  be.  shown  in 
Appendix  B,  these  terms  are  omitted. 


AFFENDIX  B 


DERIVATION  OF  C  AND  Prt  FUNCTIONS  FOR  ENTIRE  FLOW  FIELD 
V 

Relationships  for  and  Prt  are  first  developed  in 

streamline  co-ordinates  for  the  recirculating  region.  Figure  (Bl) 
shows  a  schematic  of  the  streamline  co-ordinate  system  In  the 
present  formulation  the  following  assunptions  of  Leschziner  and 
Rodi  164]  are  used. 
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Se-writing  Bgn.  (3.4.27)  for  each  ccnjnnent  in  (s,n)  co-ordinate 
system. 
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Substitution  of  Eqns.  (B6)-(Ril)  into  Bgns.  (B2)-(B4)  results  in 
the  following, 
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Bgns.  (B12)-(B14)  represent  three  algebraic  equations  with  three 
unknown  stresses,  u^u^,  u|  and  .  Elimination  of  and  u^  in 
Eqn.  (EL2)  results  in  the  expression  for  as  given  by  Bgn. 

(3.4.32).  Ocrparing  Bqn.  (3.4.32)  with  the  Boussinesq 
approximation,  Bgn.  (B5),  results  in  the  functional  form  for  , 
as  given  by  Bgn.  (3.4,36). 

Following  a  similar  procedure  as  outlined  above,  the 
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functional  representation  for  Prt  is  obtained.  Writing  Eqn. 
(3.4.31)  in  (s,n)  co-ordinates,  for  both  fluxes, 
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The  Boussinesq  approximation  is. 
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Elimination  of  u^  e‘  in  Eqn.  (B16)  by  Eqn.  (B15)  and  comparing  the 
result  with  Eqn.  (B17)  give  the  functional  form  of  Prf  TO  isolate 
the  main  contributions  of  curvature'  effects  on  Prt,  the  following 
assumption  is  made. 


9T  ,  9T 

Ss  /  sr 


0 


With  this  assumption  the  expression  for  Prt  is  obtained  as  given  by 
Eqn.  (3.4.37) .  Expressions  are  now  developed  for  Rc  and  /9n . 
Using  the  definition  of  the  radius  of  curvature, 
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Using  the  functional  relationship  for 
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Substituting 

these 

expressions 

into 

Eqn. 

(B18) 

result. 

expressions  for  ^  as  given  by  Bqn.  (3.4.12). 

3  Us/  3  n  is  determined  in  the  following  manner.  Transforming 
from  (x,y)  to  (s,n)  co-ordinates, 


S. «  S.  C..  C. 5 
k£  ij  lk  j  Z 

where  the  C’s  are  the  direction  cosines, 
in  (n,s)  system. 


(B19) 

Rewriting  equation  (BX9) 


Sns  "  IsinC2e)l  syy  '  +  Sxy  cos^  (B20) 


Talcing  S  to  be  the  shear  stresses. 
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In  the  above  expression  j“0,l  correspond  to  the  rearward  facing 
step  and  pipe  expansion  geometries,  respectively.  In  the 
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A  similar  analysis  for  Pr^  is  performed,  the  final  result  is  given 
fay  Bqn.  (3.4.39). 


APPENDIX  C 


SIMULATION  OF  INITIAL  CONDITIONS  AND  TREATMENT  OF  NEAR-WALL  FLOWS 
D.l  Simulation  of  Initial  Conditions 

Figure  (Cl)  illustrates  the  schematic  procedure  for  the 
simulation  of  approach  boundary  layer  thicknesses  for  the  Aung  and 
Goldstein  [8]  experiment.  The  experiment  was  performed  with  an 
unheated  starting  length,  The  actual  dimensions  are, 

S;  «*  15  cm.  rS'  -  x",  £  =  *1  -  X2 

xi  *  45 *.5  cm. 

X2  *  30.5  cm- 
xg  *  10  cm. 

Eton  laminar  flat  plate  boundary  layer  theory  [l03j 


Equating  the  turbulent  boundary  layer  thicknesses,  Ftps.  (3.6.1) 
and  (3.6.2)  with  Eqns  (Cl)  and  (C2)  results  in  expressions  for  the 
respective  starting  lengths.  The  effective  starting  lengths  are, 
xx  »  xs  +  xi 
*2  *  xs  +  %2 

Inserting  these  lengths  into  Eqns.  (3.6.1)  and  (3.6.2)  results  in 
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the  momentum  and  thermal  boundary  layer  thicknesses  at  the  step. 

D.2  Wall  Boundary  Condition 

Figure  (C2)  illustrates  the  mean  velocity  distribution  in  a 
boundary  layer  [12] ,  also  indicated  are  the  names  associated  with 
each  region.  A  brief  mathematical  description  is  presented  for  the 
functional  behavior  of  the  mean  velocity  in  each  of  the  regions, 
ihe  analysis  is  formulated  for  the  u-velocity;  however,  similar 
procedures  can  be  carried  out  for  the  v- velocity  and  temperature. 

1.  Linear  Sublayer 

The  viscous  sublayer  is  approximately  equal  to  o.l  percent 
of  the  total  boundary  layer  thickness,  this  corresponds  to  y+  <  5. 

^  i  - — 

n.c  twun  oii^aL  uuLcoa  10  / 


T  ♦  X 
i  t 


where. 


(laminar) 


9u 

y  ^7 


CC3) 


tt  (turbulent)  =  -  p  u'v'  (Boussinesq  Approximation) 

Ihe  munentum  equation  for  a  two-dimensional,  incompress  ible , 
turbulent  boundary  layer  with  zero  pressure  gradient  is, 


/  LINEAR  SUBLAYER 
i!  BUFFER  ZONE 
III  LOG-LAW  REGION 
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Figure  C2.  Mean  Velocity  Distribution  in  Turbulent  Bourdary 


Layer 
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u|H.  +  v  -  Ip. 

3y  p  3y 

where, 

u*v«0,  at  y  *  Q 
It  can  be  shewn  that. 


(C4) 


32t' 

ay2 


y=0 


Thus,  for  a  snail  distance  away  frem  the  wall. 


w  -° 


This  implies  that  the  region  close  to  the  wall  is  a  constant  shear 
stress  region.  Be-writing  Eqn.  (C3), 


x  du  '  ’  w  _  2 

1  ■  V  T-  -  U  V  s  —  =  u 

p  dy  p  x 


CCS) 


In  the  viscous  sublayer  region,  Bgn.  (C5)  reduces  to  (since  u'v 


I 

m 


0  at  y  «  Oj, 


lb2 


(C6) 


Integrating  Eqn.  (C6),  and  using  the  definition  of, 
y+  =  yuT  and  u+  «  u  ,  result  in, 


u+  -  y+  (C7) 

2.  Buffler  Layer 

For  this  region,  often  referred  to  as  the  transition  or  inner 
region,  y+  ranges  from  about  5  to  13.  In  this  range  both  and 
xt  are  important.  Fran  the  mixing  length  hypothesis, 


?  i 

-  p  u  v 


8u 

w 


CCS) 


Using  the  Van  Driest  [56]  length  sode  (Bgn.  (3.4.4))  and 
subsituting  into  Eqn.  (C8),  yields. 


v 


du 

Oy 


(icy)2  [l  -  exp 


For  a  (y+)  =  (icy'")2  [  1  -  exp  C-y/A*)]2 

and,  b  *=  1 

Bgn.  (C9)  reduces  to. 


(C9) 
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or, 

du+  -b  «■  (b2  +  4a) 1/2 

"  +  2a 

dy 

or, 

u*  *  /  -  ri--TC  -*/  '  «™> 

0  b  +  J  b^  +  4a(y  )] 

Egn.  (CIO)  agrees  very  well  with  the  experimental  data  of  Klebanoff 
1112J  and  Laufer  [113J  in  the  buffer  zone  [12] . 

3.  log-wall  Region 

For  this  region,  sometimes  referred  to  as  the  fully 
turbulent  part  of  the  inner  region,  y*  >13.  Egn.  (C3),  reduces 


(Cll) 


By  substituting  the  Boussinesq  approximation,  Bgn.  (3.2.1)  and  the 
expression  for  the  turbulent  viscosity,  yt,  Bgn.  (3.4.22),  into 
Bgn.  (C  11),  result  in 


<  uT  y 


du 

3y 


(C12) 


-  Integrating  Eqn.  (C  12)  results  in. 
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Figure  3.  Effect  of  Approach  Boundary  Layer  Thickness 
On  Reattachrent  length 
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Figure  4.  Geccnetrical  Representation  cf  Iv^-DimensioiAl 
Back'-Step  and  Sudden,  Pipe  Expansion 
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Ocnparisgn  of  Streanwise  Nusselt  Number  Downs  bream  of  Step  by  Quadratic 
and  Hybrid  Differencing  Schemes  for  Present  Iirproved  k-c  Model 
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PREDICTION,  PRESENT  IMPROVED  k -C  MODEL, 

PREDICTION ,  PRESENT  IMPROVED  k-€  MODEL  WITH 
SENSITIVITY  . CORRECTION  OF  LESCHZINER  AND 
ROD I  [64]  (EQU  3-4-10)  X„/«=6.0 
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SIMULATED  INITIAL  CONDITIONS,  Xn/H=5.8 
INFLUENCE  OF  TRIP ,  Xh/H=5.3 
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Experiment:  Eaton  &  Johnston  [  20  | 
Present  Prediction  (Tuo-Pass  Procedure) 


Figure  25. 
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Figure  27.  Tarperature  Distribution  Downstream  of  Step 
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.  EXPERIMENT :  MOSS  &  BAKER  QsQ 
— PRESENT  PREDICTION  (TWO -PASS  PROCEDURE) 


two-pass  ppccbdmb  mfsfwr 

IMPROVED  STANDARD  ‘A-e  MODEL  predictions 

IMPROVED  LOW -Re. No.  MODEL 


PRESENT  PREDICTION ,  TWO-PASS  PROCEDURE 


PRESENT PREDICTION,  TWO -PASS  PROCEDURE 


EXPERIMENT:  ZEMANICK  (  DOUGALL  [39] 

CHSENG  f  LAUNDER  [59] 

PRESENT  PREDICTION ,  71 WJ- PASS  PROCEDURE 


EXPERIMENT:  ZEMANICK  (  DOUGAUL  [39] 
CHENG  $  LAUNDER  [59] 

PRESENT  PREDICTION *  TWO -PASS  PROCEDURE 


Nu. 


max 


/ 


/ 


< 


/ 


/o 

1U 


0.54 


0  ZEMANICK  i  DOUG  ALL  USD 
- PRESENT  PREDICTION,  A/ifo0~  Re0?5 

10  i  1 _ I L_l _ I - ! - 1 — L_l 

IOl  .  10 4  /os 

Be 


Figure  38.  Variation  of  Maximum  Nusselt  Number  with 

Reynolds  Number  Downstream  of  a  Sudden 
Pipe  Expansion 
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Figure  40  o  Cross-Stream  Variation  of  Nusselt  Number 

on  the  Downstream  Face  of  the  Pipe  Expansion 
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Figure  41.  Croi^-Str&cr.  Variation  of  Kosaelt  Nlsnber  on 
the  Dowrastnasm  Faoe  of  the  Pipe  Expansion 


©  EXPERIMENT:  SPARROW  $  O'BRIEN  1 4Q 
_ PRESENT  PREDICTION 


(2y-<!)/(D-d) 


Figure  42.  Cross-Stream  Variation  of  Nusselt  Number  on  tha 
Downstream  Face  of  the  Pipe  Expansion 
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Figure  43.  Variation  of  Average  Nussslt  Number  with  Reynolds  Number 
for  the  Dwnstream  Face  of  a  Sudden  Pipe  ESqpansion 
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Figure  48.  Streamwi.se  Nusselt  Nurrber  for  Back-Step  Normalized 
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